Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,521,284 A→G D486G (GAT→GGT)  igaA → intracellular growth attenuator protein IgaA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,521,2840AG97.3% 261.9 / ‑1.6 73D486G (GAT→GGT) igaAintracellular growth attenuator protein IgaA
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (32/39);  total (33/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.17e-01

AACTGGTCAACAAAGCGACGGCATTGACTGAAGCAGTTAATCGCCAGCTGCACCCTAAACCGGAAGATGAATCTCGCGTCAGTGCCTCATTACGTTCAGCAATTCAAAAATCCGGCATGGTATTGCTTGATGATTTTGGCGACATTGTACTGAAGACAGCGGATTTATGTTCTGCCAAAGATGACTGTGTGCGACTGAAAAATGCGCTGGTCAATCTCGGCAACAGTAAAGACTGGGACGCGCTGGTAAAACGCGCCAACGCCGGG  >  NZ_CP009273/3521152‑3521417
                                                                                                                                    |                                                                                                                                     
aaCTGGTCAACAAAGCGACGGCATTGACTGAAGCAGTTAATCGCCAGCTGCACCCTAAACCGGAAGATGAATCTCGCGTCAGTGCCTCATTACGTTCAGCAATTCAAAAATCCGGCATGGTATTGCTTGATGGttt                                                                                                                                    <  1:765442/136‑1 (MQ=255)
aaCTGGTCAACAAAGCGACGGCATTGACTGAAGCAGTTAATCGCCAGCTGCACCCTAAACCGGAAGATGAATCTCGCGTCAGTGCCTCATTACGTTCAGCAATTCAAAAATCCGGCATGGTATTGCTTGATGGTTTTgg                                                                                                                                 >  2:576377/1‑139 (MQ=255)
    ggTCAACAAAGCGACGGCATTGACTGAAGCAGTTAATCGCCAGCTGCACCCTAAACCGGAAGATGAATCTCGCGTCAGTGCCTCATTACGTTCAGCAATTCAAAAATCCGGCATGGTATTGCTTGATGGTTTTGGCGa                                                                                                                              >  1:655630/1‑138 (MQ=255)
     gTCAACAAAGCGACGGCATTGACTGAAGCAGTTAATCGCCAGCTGCACCCTAAACCGGAAGATGAATCTCGCGTCAGTGCCTCATTACGTTCAGCAATTCAAAAATCCGGCATGGTATTGCTTGATGGTTTTGGCGACa                                                                                                                            <  2:281948/139‑1 (MQ=255)
     gTCAACAAAGCGACGGCATTGACTGAAGCAGTTAATCGCCAGCTGCACCCTAAACCGGAAGATGAATCTCGCGTCAGTGCCTCATTACGTTCAGCAATTCAAAAATCCGGCATGGTATTGCTTGATGGTTTTGGCGACa                                                                                                                            <  1:126056/139‑1 (MQ=255)
     gTCAACAAAGCGACGGCATTGACTGAAGCAGTTAATCGCCAGCTGCACCCTAAACCGGAAGATGAATCTCGCGTCAGTGCCTCATTACGTTCAGCAATTCAAAAATCCGGCATGGTATTGCTTGATGGTTTTGGCGACa                                                                                                                            <  2:348054/139‑1 (MQ=255)
       caacaaAGCGACGGCATTGACTGAAGCAGTTAATCGCCAGCTGCACCCTAAACCGGAAGATGAATCTCGCGTCAGTGCCTCATTACGTTCAGCAATTCAAAAATCCGGCATGGTATTGCTTGATGGTTTTGGCGACAtt                                                                                                                          >  2:768698/1‑139 (MQ=255)
             aGCGACGGCATTGACTGAAGCAGTTAATCGCCAGCTGCACCCTAAACCGGAAGATGAATCTCGCGTCAGTGCCTCATTACGTTCAGCAATTCAAAAATCCGGCATGGTATTGCTTGATGGTTTTGGCGACATTGTACt                                                                                                                     >  1:1138282/1‑138 (MQ=255)
               cGACGGCATTGACTGAAGCAGTTAATCGCCAGCTGCACCCTAAACCGGAAGATGAATCTCGCGTCAGTGCCTCATTACGTTCAGCAATTCAAAAATCCGGCATGGTATTGCTTGATGGTTTTGGCGACATTGTACTg                                                                                                                    >  1:1254541/1‑137 (MQ=255)
               cGACGGCATTGACTGAAGCAGTTAATCGCCAGCTGCACCCTAAACCGGAAGATGAATCTCGCGTCAGTGCCTCATTACGTTCAGCAATTCAAAAATCCGGCATGGTATTGCTTGATGGTTTTGGCGACATTGTACTg                                                                                                                    >  2:546960/1‑137 (MQ=255)
               cGACGGCATTGACTGAAGCAGTTAATCGCCAGCTGCACCCTAAACCGGAAGATGAATCTCGCGTCAGTGCCTCATTACGTTCAGCAATTCAAAAATCCGGCATGGTATTGCCTGATGGTTTTGGCGACATTGTACTg                                                                                                                    >  2:875672/1‑137 (MQ=255)
                   ggCATTGACTGAAGCAGTTAATCGCCAGCTGCACCCTAAACCGGAAGATGAATCTCGCGTCAGTGCCTCATTACGTTCAGCAATTCAAAAATCCGGCATGGTATTGCTTGATGGTTTTGGCGACATTGTACTGAAGACa                                                                                                              <  1:247552/139‑1 (MQ=255)
                    gCATTGACTGAAGCAGTTAATCGCCAGCTGCACCCTAAACCGGAAGATGAATCTCGCGTCAGTGCCTCATTACGTTCAGCAATTCAAAAATCCGGCATGGTATTGCTTGATGGTTTTGGCGACATTGTACTGAAGACAg                                                                                                             <  1:80943/139‑1 (MQ=255)
                    gCATTGACTGAAGCAGTTAATCGCCAGCTGCACCCTAAACCGGAAGATGAATCTCGCGTCAGTGCCTCATTACGTTCAGCAATTCAAAAATCCGGCATGGTATTGCCTGATGGTTTTGGCGACATTGTACTGAAGACAg                                                                                                             <  1:875672/139‑1 (MQ=255)
                        tGACTGAAGCAGTTAATCGCCAGCTGCACCCTAAACCGGAAGATGAATCTCGCGTCAGTGCCTCATTACGTTCAGCAATTCAAAAATCCGGCATGGTATTGCTTGATGGTTTTGGCGACATTGTACTGAAGACAGCGGa                                                                                                         >  1:1257174/1‑139 (MQ=255)
                               aGCAGTTAATCGCCAGCTGCACCCTAAACCGGAAGATGAATCTCGCGTCAGTGCCTCATTACGTTCAGCAATTCAAAAATCCGGCATGGTATTGCTTGATGGTTTTGGCg                                                                                                                               <  2:243862/110‑1 (MQ=255)
                               aGCAGTTAATCGCCAGCTGCACCCTAAACCGGAAGATGAATCTCGCGTCAGTGCCTCATTACGTTCAGCAATTCAAAAATCCGGCATGGTATTGCTTGATGGTTTTGGCg                                                                                                                               >  1:243862/1‑110 (MQ=255)
                               aGCAGTTAATCGCCAGCTGCACCCTAAACCGGAAGATGAATCTCGCGTCAGTGCCTCATTACGTTCAGCAATTCAAAAATCCGGCATGGTATTGCTTGATGGTTTTGGCGACATTGTACTGAAGACAGCGGATTTATGt                                                                                                  >  1:401806/1‑139 (MQ=255)
                               aGCAGTTAATCGCCAGCTGCACCCTAAACCGGAAGATGAATCTCGCGTCAGTGCCTCATTACGTTCAGCAATTCAAAAATCCGGCATGGTATTGCTTGATGGTTTTGGCGACATTGTACTGAAGACAGCGGATTTATGt                                                                                                  >  2:871161/1‑139 (MQ=255)
                                 cAGTTAATCGCCAGCTGCACCCTAAACCGGAAGATGAATCTCGCGTCAGTGCCTCATTACGTTCAGCAATTCAAAAATCCGGCATGGTATTGCTTGATGGTTTTGGCGACATTGTACTGAAGACAGCGGATTTATGTTc                                                                                                <  1:871161/139‑1 (MQ=255)
                                   gTTAATCGCCAGCTGCACCCTAAACCGGAAGATGAATCTCGCGTCAGTGCCTCATTACGTTCAGCAATTCAAAAATCCGGCATGGTATTGCTTGATGGTTTTGGCGACATTGTACTGAAGACAGCGGATTTATGTTCTg                                                                                              <  1:392695/139‑1 (MQ=255)
                                   gTTAATCGCCAGATGCACCCTAAACCGGAAGATGAATCTCGCGTCAGTGCCTCATTACGTTCAGCAATTCAAAAATCCGGCATGGTATTGCTTGATGGTTTTGGCGACATTGTACTGAAGACAGCGGATTTATGTTCTg                                                                                              <  2:361738/139‑1 (MQ=255)
                                    ttAATCGCCAGCTGCACCCTAAACCGGAAGATGAATCTCGCGTCAGTGCCTCATTACGTTCAGCAATTCAAAAATCCGGCATGGTATTGCTTGATGGTTTTGGCGACATTGTACTGAAGACAGCGGATTTATGTTCTGc                                                                                             >  1:518953/1‑139 (MQ=255)
                                      aaTCGCCAGCTGCACCCTAAACCGGAAGATGAATCTCGCGTCAGTGCCTCATTACGTTCAGCAATTCAAAAATCCGGCATGGTATTGCTTGATGGTTTTGGCGACATTGTACTGAAGACAGCGGATTTATGTTCTGcc                                                                                            <  1:768698/138‑1 (MQ=255)
                                        tCGCCAGCTGCACCCTAAACCGGAAGATGAATCTCGCGTCAGTGCCTCATTACGTTCAGCAATTCAAAAATCCGGCATGGTATTGCTTGATGGTTTTGGCGACATTGTACTGAAGACAGCGGATTTATGTTCTGCCa                                                                                           >  2:976964/1‑137 (MQ=255)
                                          gCCAGCTGCACCCTAAACCGGAAGATGAATCTCGCGTCAGTGCCTCATTACGTTCAGCAATTCAAAAATCCGGCATGGTATTGCTTGATGGTTTTGGCGACATTGTACTGAAGACAGCGGATTTATGTTCTGCCAAAGa                                                                                       <  2:1254541/139‑1 (MQ=255)
                                             aGCTGCACCCTAAACCGGAAGATGAATCTCGCGTCAGTGCCTCATTACGTTCAGCAATTCAAAAATCCGGCATGGTATTGCTTGATGGTTTTGGCGACATTGTACTGAAGACAGCGGATTTATGTTCTGCCAAAGATGa                                                                                    <  2:1257174/139‑1 (MQ=255)
                                              gCTGCACCCTAAACCGGAAGATGAATCTCGCGTCAGTGCCTCATTACGTTCAGCAATTCAAAAATCCGGCATGGTATTGCTTGATGGTTTTGGCGACATTGTACt                                                                                                                     >  2:51024/1‑105 (MQ=255)
                                              gCTGCACCCTAAACCGGAAGATGAATCTCGCGTCAGTGCCTCATTACGTTCAGCAATTCAAAAATCCGGCATGGTATTGCTTGATGGTTTTGGCGACATTGTACt                                                                                                                     <  1:51024/105‑1 (MQ=255)
                                              gCTGCACCCTAAACCGGAAGATGAATCTCGCGTCAGTGCCTCATTACGTTCAGCAATTCAAAAATCCGGCATGGTATTGCTTGATGGTTTTGGCGACATTGTACTGAAGACAGCGGATTTATGTTCTGCCAAAGCTGAc                                                                                   >  2:342671/1‑139 (MQ=255)
                                              gCTGCACCCTAAACCGGAAGATGAATCTCGCGTCAGTGCCTCATTACGTTCAGCAATTCAAAAATCCGGCATGGTATTGCTTGATGGTTTTGGCGACATTGTACTGAAGACAGCGGATTTATGTTCTGCCAAAGATGAc                                                                                   >  2:342688/1‑139 (MQ=255)
                                                    cccTAAACCGGAAGATGAATCTCGCGTCAGTGCCTCATTACGTTCAGCAATTCAAAAATCCGGCATGGTATTGCTTGATGGTTTTGGCGACATTGTACTGAAGACAGCGGATTTATGTTCTGCCAAAGATGACtgtgtg                                                                             <  1:976964/139‑1 (MQ=255)
                                                    cccTAAACCGGAAGATGAATCTCGCGTCAGTGCCTCATTACGTTCAGCAATTCAAAAATCCGGCATGGTATTGCTTGATGGTTTTGGCGACATTGTACTGAAGACAGCGGATTTATGTTCTGCCAAAGATGACtgtgtg                                                                             <  2:132366/139‑1 (MQ=255)
                                                       tAAACCGGAAGATGAATCTCGCGTCAGTGCCTCATTACGTTCAGCAATTCAAAAATCCGGCATGGTATTGCTTGATGGTTTTGGCGACa                                                                                                                            <  2:286805/89‑1 (MQ=255)
                                                       tAAACCGGAAGATGAATCTCGCGTCAGTGCCTCATTACGTTCAGCAATTCAAAAATCCGGCATGGTATTGCTTGATGGTTTTGGCGACa                                                                                                                            >  1:286805/1‑89 (MQ=255)
                                                         aaCCGGAAGATGAATCTCGCGTCAGTGCCTCATTACGTTCAGCAATTCAAAAATCCGGCATGGTATTGCTTGATGGTTTTGGCGACATTGTACTGAAGACAGCGGATTTATGTTCTGCCAAAGATGACTGTGTGCGa                                                                          <  2:270308/137‑1 (MQ=255)
                                                              gAAGATGAATCTCGCGTCAGTGCCTCATTACGTTCAGCAATTCAAAAATCCGGCATGGTATTGCTTGATGGTTTTGGCGACATTGTACTGAAGACAGCGGATTTATGTTCTGCCAAAGATGACTGTGTGCGACTGaaa                                                                    >  2:1064281/1‑138 (MQ=255)
                                                               aaGATGAATCTCGCGTCAGTGCCTCATTACGTTCAGCAATTCAAAAATCCGGCATGGTATTGCTTGATGGTTTTgg                                                                                                                                 <  2:969657/76‑1 (MQ=255)
                                                               aaGATGAATCTCGCGTCAGTGCCTCATTACGTTCAGCAATTCAAAAATCCGGCATGGTATTGCTTGATGGTTTTgg                                                                                                                                 >  1:969657/1‑76 (MQ=255)
                                                               aaGATGAATCTCGCGTCAGTGCCTCATTACGTTCAGCAATTCAAAAATCCGGCATGGTATTGCTTGATGGTTTTGGCGACATTGTACTGAAGACAGCGGATTTATGTTCTGCCAAAGATGACTGTGTGCGACTGaaaaa                                                                  >  1:716947/1‑139 (MQ=255)
                                                               aaGATGAATCTCGCGTCAGTGCCTCATTACGTTCAGCAATTCAAAAATCCGGCATGGTATTGCTTGATGGTTTTGGCGACATTGTACTGAAGACAGCGGATTTATGTTCTGCCAAAGATGACTGTGTGCGACTGaaaaa                                                                  >  2:575739/1‑139 (MQ=255)
                                                                   tGAATCTCGCGTCAGTGCCTCATTACGTTCAGCAATTCAAAAATCCGGCATGGTATTGCTTGATGGTTTTGGCGACATTGTACTGAAGACAGCGGATTTATGTTCTGCCAAAGATGACTGTGTGCGACTGAAAAATgcg                                                              <  1:880486/139‑1 (MQ=255)
                                                                      aTCTCGCGTCAGTGCCTCATTACGTTCAGCAATTCAAAAATCCGGCATGGTATTGCTTGATGGTTTTGGCGACATTGTACTGAAGACAGCGGATTTATGTTCTGCCAAAGATGACTGTGTGCGACTGAAAAATGCGCTg                                                           >  2:171627/1‑139 (MQ=255)
                                                                         tcGCGTCAGTGCCTCATTACGTTCAGCAATTCAAAAATCCGGCATGGTATTGCTTGATGGTTTTGGCGACATTGTACTGAAGACAGCGGATTTATGTTCTGCCAAAGATGACTGTGTGCGACTGAAAAATGCGCTGGt                                                         >  1:516227/1‑138 (MQ=255)
                                                                           gcgTCAGTGCCTCATTACGTTCAGCAATTCAAAAATCCGGCATGGTATTGCTTGATGGTTTTGGCGACATTGTACTGAAGACAGCGGATTTATGTTCTGCCAAAGATGACTGTGTGCGACTGAAAAATGCGCTGGt                                                         <  2:1032605/136‑1 (MQ=255)
                                                                           gcgTCAGTGCCTCATTACGTTCAGCAATTCAAAAATCCGGCATGGTATTGCTTGATGGTTTTGGCGACATTGTACTGAAGACAGCGGATTTATGTTCTGCCAAAGATGACTGTGTGCGACTGAAAAATGCGCTGGt                                                         <  1:577886/136‑1 (MQ=255)
                                                                           gcgTCAGTGCCTCATTACGTTCAGCAATTCAAAAATCCGGCATGGTATTGCTTGATGGTTTTGGCGACATTGTACTGAAGACAGCGGATTTATGTTCTGCCAAAGATGACTGTGTGCGACTGAAAAATGCGCTGGt                                                         <  1:468620/136‑1 (MQ=255)
                                                                           gcgTCAGTGCCTCATTACGTTCAGCAATTCAAAAATCCGGCATGGTATTGCTTGATGGTTTTGGCGACATTGTACTGAAGACAGCGGATTTATGTTCTGCCAAAGATGACTGTGTGCGACTGAAAAATGCGCTGGTCa                                                       <  1:1094664/138‑1 (MQ=255)
                                                                            cgTCAGTGCCTCATTACGTTCAGCAATTCAAAAATCCGGCATGGTATTGCTTGATGGTTTTGGCGACATTGTACTGAAGACa                                                                                                              <  1:1063785/82‑1 (MQ=255)
                                                                            cgTCAGTGCCTCATTACGTTCAGCAATTCAAAAATCCGGCATGGTATTGCTTGATGGTTTTGGCGACATTGTACTGAAGACa                                                                                                              >  2:1063785/1‑82 (MQ=255)
                                                                                 gTGCCTCATTACGTTCAGCAATTCAAAAATCCGGCATGGTATTGCTTGATGGTTTTgg                                                                                                                                 >  1:411567/1‑58 (MQ=255)
                                                                                 gTGCCTCATTACGTTCAGCAATTCAAAAATCCGGCATGGTATTGCTTGATGGTTTTgg                                                                                                                                 <  2:411567/58‑1 (MQ=255)
                                                                                 gTGCCTCATTACGTTCAGCAATTCAAAAATCCGGCATGGTATTGCTTGATGGTTTTGGCGACATTGTACTGAAGACAGCGGATTTATGTTCTGCCAAAGATGACTGTGTGCGACTGAAAAATGCGCTGGTCAATCTCgg                                                <  2:401806/139‑1 (MQ=255)
                                                                                  tGCCTCATTACGTTCAGCAATTCAAAAATCCGGCATGGTATTGCTTGATGGTTTTGGCGACATTGTACTGAAGACAGCGGATTTATGTTCTGCCAAAGATGACTGTGTGCGACTGAAAAATGCGCTGGTCAATCTCGGc                                               <  2:518953/139‑1 (MQ=255)
                                                                                    ccTCATTACGTTCAGCAATTCAAAAATCCGGCATGGTATTGCTTGATGGTTTTGGCGACATTGTACTGAAGACAGCGGATTTATGTTCTGCCAAAGATGACTGTGTGCGACTGAAAAATGCGCTGGTCAATCTCGGCa                                              <  2:435387/138‑1 (MQ=255)
                                                                                     cTCATTACGTTCAGCAATTCAAAAATCCGGCATGGTATTGCTTGATGGTTTTGGCGACATTGTACTGAAGACAGCGGATTTATGTTCTGCCAAAGATGACTGTGTGCGACTGAAAAATGCGCTGGTCAATCTCGGCaa                                             >  1:739082/1‑138 (MQ=255)
                                                                                        aTTACGTTCAGCAATTCAAAAATCCGGCATGGTATTGCTTGATGGTTTTGGCGACATTGTACTGAAGACAGCGGATTTATGTTCTGCCAAAGATGACTGTGTGCGACTGAAAAATGCGCTGGTCAATCTCGGCAACAGt                                         <  2:507213/139‑1 (MQ=255)
                                                                                            cGTTCAGCAATTCAAAAATCCGGCATGGTATTGCTTGATGATTTTGGCGACATTGTACTGAAGACAGCg                                                                                                           <  2:887084/69‑1 (MQ=255)
                                                                                            cGTTCAGCAATTCAAAAATCCGGCATGGTATTGCTTGATGATTTTGGCGACATTGTACTGAAGACAGCg                                                                                                           >  1:887084/1‑69 (MQ=255)
                                                                                                cAGCAATTCAAAAATCCGGCATGGTATTGCTTGATGGTTTTGGCGACATTGTACTGAAGACAGCGGATTTATGTTCTGCCAAAGATGACTGTGTGCGACTGAAAAATGCGCTGGTCAATCTCGGc                                               <  2:1290542/125‑1 (MQ=255)
                                                                                                cAGCAATTCAAAAATCCGGCATGGTATTGCTTGATGGTTTTGGCGACATTGTACTGAAGACAGCGGATTTATGTTCTGCCAAAGATGACTGTGTGCGACTGAAAAATGCGCTGGTCAATCTCGGc                                               >  1:1290542/1‑125 (MQ=255)
                                                                                                     aTTCAAAAATCCGGCATGGTATTGCTTGATGGTTTTGGCGACATTGTACTGAAGACAGCGGATTTATGTTCTGCCAAAGATGACTGTGTGCGACTGAAAAATGCGCTGGTCAATCTCGGCAACAGTAAAGACTGGGAcg                            <  1:576377/139‑1 (MQ=255)
                                                                                                        cAAAAATCCGGCATGGTATTGCTTGATGGTTTTGGCGACATTGTACTGAAGACAGCGGATTTATGTTCTGCCAAAGATGACTGTGTGCGACTGAAAAATGCGCTGGTCAATCTCGGCAACAGt                                         <  1:173824/123‑1 (MQ=255)
                                                                                                        cAAAAATCCGGCATGGTATTGCTTGATGGTTTTGGCGACATTGTACTGAAGACAGCGGATTTATGTTCTGCCAAAGATGACTGTGTGCGACTGAAAAATGCGCTGGTCAATCTCGGCAACAGt                                         >  2:173824/1‑123 (MQ=255)
                                                                                                        cAAAAATCCGGCATGGTATTGCTTGATGGTTTTGGCGACATTGTACTGAAGACAGCGGATTTATGTTCTGCCAAAGATGACTGTGTGCGACTGAAAAATGCGCTGGTCAATCTCGGCAACAGTAAAGACTGGGa                              <  1:793934/134‑1 (MQ=255)
                                                                                                        cAAAAATCCGGCATGGTATTGCTTGATGGTTTTGGCGACATTGTACTGAAGACAGCGGATTTATGTTCTGCCAAAGATGACTGTGTGCGACTGAAAAATGCGCTGGTCAATCTCGGCAACAGTAAAGACTGGGa                              >  2:793934/1‑134 (MQ=255)
                                                                                                        cAAAAATCCGGCATGGTATTGCTTGATGGTTTTGGCGACATTGTACTGAAGACAGCGGATTTATGTTCTGCCAAAGATGACTGTGTGCGACTGAAAAATGCGCTGGTCAATCTCGGCAACAGTAAAGACTGGGAcgcg                          <  1:752440/138‑1 (MQ=255)
                                                                                                           aaaTCCGGCATGGTATTGCTTGATGGTTTTGGCGACATTGTACTGAAGACAGCGGATTTATGTTCTGCCAAAGATGACTGTGTGCGACTGAAAAATGCGCTGGTCAATCTCGGCAACAGTAAAGACTGGGAcgcg                          <  2:469983/135‑1 (MQ=255)
                                                                                                           aaaTCCGGCATGGTATTGCTTGATGGTTTTGGCGACATTGTACTGAAGACAGCGGATTTATGTTCTGCCAAAGATGACTGTGTGCGACTGAAAAATGCGCTGGTCAATCTCGGCAACAGTAAAGACTGGGAcgcg                          >  1:469983/1‑135 (MQ=255)
                                                                                                                             cTTGATGGTTTTGGCGACATTGTACTGAAGACAGCGGATTTATGTTCTGCCAAAGAt                                                                                      >  1:954871/1‑57 (MQ=255)
                                                                                                                             cTTGATGGTTTTGGCGACATTGTACTGAAGACAGCGGATTTATGTTCTGCCAAAGAt                                                                                      <  2:954871/57‑1 (MQ=255)
                                                                                                                             cTTGATGGTTTTGGCGACATTGTACTGAAGACAGCGGATTTATGTTCTGCCAAAGATGACTGTGTGCGACTGAAAAATGCGCTGGTCAATCTCGGCAACAGTAAAGACTGGGACGCGCTGGTaaaa                 >  2:998633/1‑126 (MQ=255)
                                                                                                                             cTTGATGGTTTTGGCGACATTGTACTGAAGACAGCGGATTTATGTTCTGCCAAAGATGACTGTGTGCGACTGAAAAATGCGCTGGTCAATCTCGGCAACAGTAAAGACTGGGACGCGCTGGTaaaa                 <  1:998633/126‑1 (MQ=255)
                                                                                                                                gatggtTTTGGCGACATTGTACTGAAGACAGCGGATTTATGTTCTGCCAAAGATGACTGTGTGCGACTGAAAAATGCGCTGGTCAATCTCGGCAACAGTAAAGACTGGGACGCGCTGGTAAAACGCGCCAACGCCggg  >  1:76203/1‑138 (MQ=255)
                                                                                                                                gatggtTTTGGCGACATTGTACTGAAGACAGCAGATTTATGTTCTGCCAAAGATGACTGTGTGCGACTGAAAAATGCGCTGGTCAATCTCGGCAACAGTAAAGACTGGGACGCGCTGGTAAAACGCGCCAACGCCggg  <  2:1138282/138‑1 (MQ=255)
                                                                                                                                 atggtTTTGGCGACATTGTACTGAAGACAGCGGATTTATGTTCTGCCAAAGATGACTGTGTGCGACTGAAAAATGCGCTGGTCAATCTCGGCAACAGTAAAGACTGGGa                              <  2:401245/109‑1 (MQ=255)
                                                                                                                                 atggtTTTGGCGACATTGTACTGAAGACAGCGGATTTATGTTCTGCCAAAGATGACTGTGTGCGACTGAAAAATGCGCTGGTCAATCTCGGCAACAGTAAAGACTGGGa                              >  1:401245/1‑109 (MQ=255)
                                                                                                                                 atggtTTTGGCGACATTGTACTGAAGACAGCGGATTTATGTTCTGCCAAAGATGACTGTGTGCGACTGAAAAATGCGCTGGTCAATCTCGGCAACAGTAAAGACTGGGACGCGCTGGTAAAACGCGCCAACGCCggg  >  1:1052817/1‑137 (MQ=255)
                                                                                                                                    |                                                                                                                                     
AACTGGTCAACAAAGCGACGGCATTGACTGAAGCAGTTAATCGCCAGCTGCACCCTAAACCGGAAGATGAATCTCGCGTCAGTGCCTCATTACGTTCAGCAATTCAAAAATCCGGCATGGTATTGCTTGATGATTTTGGCGACATTGTACTGAAGACAGCGGATTTATGTTCTGCCAAAGATGACTGTGTGCGACTGAAAAATGCGCTGGTCAATCTCGGCAACAGTAAAGACTGGGACGCGCTGGTAAAACGCGCCAACGCCGGG  >  NZ_CP009273/3521152‑3521417

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: