Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,530,069 T→A 13.2% intergenic (+31/‑137) PP_4018 → / → topB acetyltransferase/DNA topoisomerase III

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,530,0690TA13.2% 252.2 / 28.0 93intergenic (+31/‑137)PP_4018/topBacetyltransferase/DNA topoisomerase III
Reads supporting (aligned to +/- strand):  ref base T (43/37);  new base A (5/7);  total (49/44)
Fisher's exact test for biased strand distribution p-value = 5.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.53e-01

TCAGCGCCCAGGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATGCGCCTGT  >  NC_002947/4529920‑4530215
                                                                                                                                                     |                                                                                                                                                  
tCAGCGCCCAGGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGtt                                                                                                                                                    <  1:1393223/150‑1 (MQ=255)
   gcgcCCAGGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTAc                                                                                                                                                 <  2:1270839/150‑1 (MQ=255)
    cgcCCAGGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTAcc                                                                                                                                                <  2:2052023/150‑1 (MQ=255)
     gcCCAGGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTAccc                                                                                                                                               <  2:120172/150‑1 (MQ=255)
      cccAGGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTAccag                                                                                                                                              >  2:1421602/1‑148 (MQ=255)
       ccAGGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCcgc                                                                                                                                             >  1:1623439/1‑150 (MQ=255)
           gTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCGGCAGGGCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAag                                                                                                                                         >  1:154530/1‑150 (MQ=255)
           gTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTAAACCCGCGAag                                                                                                                                         <  1:162904/150‑1 (MQ=255)
            tGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTAAACCCTCGAaga                                                                                                                                        <  1:795318/150‑1 (MQ=255)
             gCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAagag                                                                                                                                       >  2:19110/1‑150 (MQ=255)
             gCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTAAACCCGCGAagag                                                                                                                                       <  1:193165/150‑1 (MQ=255)
               acgccacgccGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCGCCTGAGCCGGCCTCTTCGCGGGTAAACCCGCGaa                                                                                                                                          >  1:916056/1‑145 (MQ=38)
               acgccacgccGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCGCCTGAGCCGGCCTCTTCGCGGGTAAACCCGCGaa                                                                                                                                          <  2:916056/145‑1 (MQ=38)
               acgccacgccGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTAAACCCGCGaa                                                                                                                                          >  1:5642/1‑145 (MQ=255)
               acgccacgccGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTAAACCCGCGaa                                                                                                                                          <  2:5642/145‑1 (MQ=255)
                   cacgccGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGc                                                                                                                                 >  1:742314/1‑150 (MQ=255)
                       ccGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAg                                                                                                                             >  2:987667/1‑150 (MQ=255)
                        cGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAgg                                                                                                                            <  1:1421602/150‑1 (MQ=255)
                               tGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGAtata                                                                                                                     >  2:181071/1‑150 (MQ=255)
                               tGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGAtata                                                                                                                     >  1:648023/1‑150 (MQ=255)
                               tGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGAtata                                                                                                                     >  1:2256708/1‑150 (MQ=255)
                                gAGCGGCTTGGCGTTCGCGTAGTCCGCGAGGAATTCCTCAAAGGTGTCACCCCGCATGTCAACATGGTGCGCCACTCGCTCGCCCGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCACGAAGAGGCCGCCTCAGGCGATATAc                                                                                                                    <  2:403243/150‑1 (MQ=255)
                                   gggCTGGGCTTTCGCGGAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCNTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATc                                                                                                                 <  1:731325/149‑1 (MQ=255)
                                       tGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCg                                                                                                             <  1:468674/150‑1 (MQ=255)
                                           cTTTCGCGTAGTCAGCGAGAAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGa                                                                                                         <  2:752689/150‑1 (MQ=255)
                                                 cgTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACaaaaa                                                                                                   <  1:181071/150‑1 (MQ=255)
                                                   tAGTCAGCGAGGAATTCCTCTAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGTTGACAAAAAAc                                                                                                 <  1:2082486/150‑1 (MQ=255)
                                                   tAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGTCTTGCGCGCCTGATCCATCCCTGGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAAc                                                                                                 <  1:1106429/150‑1 (MQ=255)
                                                   tAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAaga                                                                                                                                        <  2:515878/111‑1 (MQ=255)
                                                   tAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAaga                                                                                                                                        >  1:515878/1‑111 (MQ=255)
                                                     gTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAAccc                                                                                               <  2:2049705/150‑2 (MQ=255)
                                                     gTCAGCGAGGAATTCCTCGAAGCTGGCATCACGCAGGTCGTCATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCt                                                                                               <  1:1102346/150‑1 (MQ=255)
                                                     cccaccaaGGACTTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCCCGCCTCGCGCGCCTTATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCt                                                                                               <  1:1170072/143‑1 (MQ=255)
                                                       cAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGt                                                                                             <  1:173949/150‑1 (MQ=255)
                                                        aGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGata                                                                                                                       <  1:353049/123‑1 (MQ=255)
                                                        aGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGata                                                                                                                       >  2:353049/1‑123 (MQ=255)
                                                        aGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCGCAGGCGATATACATCCCCGCTGACAAAAAACCTGTa                                                                                            >  2:427873/1‑150 (MQ=255)
                                                         gCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTAc                                                                                           >  2:1207153/1‑150 (MQ=255)
                                                          cGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTAca                                                                                          >  2:696109/1‑150 (MQ=255)
                                                          cGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTAca                                                                                          <  1:427873/150‑1 (MQ=255)
                                                             ggAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACATGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCACTGACTAAAAAGCTTTAcatat                                                                                       >  1:98589/1‑148 (MQ=255)
                                                                atacatCGAAGCTAGAATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCAGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACAACAAAAAAAAAAAAAAACCTGTAcatccatc                                                                                    <  2:154530/145‑1 (MQ=255)
                                                                aTTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGatat                                                                                                                      <  2:1227103/116‑1 (MQ=255)
                                                                aTTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGatat                                                                                                                      >  1:1227103/1‑116 (MQ=255)
                                                                aTTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTAcatccatc                                                                                    <  1:1207153/150‑1 (MQ=255)
                                                                 ttCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACatccatcc                                                                                   >  2:1108508/1‑150 (MQ=255)
                                                                 ttCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTAAACGCGACGCCCTCTTCCTATGCGTGATAAAGTACCGCTGGCAAAAAACCTGTACatccatcc                                                                                   >  2:60514/1‑150 (MQ=255)
                                                                 ttCCTCGAAGCTGGCATCCCGCAGGTCGACACGGTGCGCGACTCGCGCGCCCGATCCATCACCTGAACTGGCCTCTTCGCGGGCCAAACCGAGAGGAGGACGGCCCAGGCGATATACAGCCCCGCGGACAAAAAACCGGTGcatacaggc                                                                                   >  2:660546/1‑147 (MQ=255)
                                                                    cTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGata                                                                                                                       <  1:1569311/111‑1 (MQ=255)
                                                                    cTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGata                                                                                                                       >  2:1569311/1‑111 (MQ=255)
                                                                     tCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCGCCTGAGCCGGCCTCTTCGCGGGTAAACCcgc                                                                                                                                             >  1:1558190/1‑88 (MQ=21)
                                                                     tCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCGCCTGAGCCGGCCTCTTCGCGGGTAAACCcgc                                                                                                                                             <  2:1558190/88‑1 (MQ=21)
                                                                     tCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAaga                                                                                                                                        <  2:502781/93‑1 (MQ=255)
                                                                     tCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAaga                                                                                                                                        >  1:502781/1‑93 (MQ=255)
                                                                     tCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGAAATATATCCCCGCTGACAAAAAACCTGTACATCCATCCAGAt                                                                               >  1:1777306/1‑150 (MQ=255)
                                                                       gAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATaa                                                                             >  1:1793853/1‑150 (MQ=255)
                                                                              gCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGc                                                                      <  2:742314/150‑1 (MQ=255)
                                                                              gCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTAAACCCGCGAaga                                                                                                                                        <  2:1882788/84‑1 (MQ=255)
                                                                              gCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTAAACCCGCGAaga                                                                                                                                        >  1:1882788/1‑84 (MQ=255)
                                                                                   ccGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACCTCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTg                                                                 <  2:1623439/150‑1 (MQ=255)
                                                                                    cGCATGTCGACAGGGGGCGCGACTCGCGGGGCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTgc                                                                <  2:1078688/150‑1 (MQ=255)
                                                                                     gCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGACGCCGGCTCAGGCGATATACATCCCCGCCGACAACAAACCTGTACATCCATCCACATAAAACTTGCCTCTgcc                                                               >  1:252511/1‑149 (MQ=255)
                                                                                             aCATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCATCTGCGCCCCCGCt                                                       >  2:2263862/1‑150 (MQ=255)
                                                                                                 ggTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTgcg                                                   >  2:1612193/1‑150 (MQ=255)
                                                                                                   tGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCa                                                 >  1:1599912/1‑150 (MQ=255)
                                                                                                    gcgcGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTAcatccatc                                                                                    <  2:404054/114‑1 (MQ=255)
                                                                                                    gcgcGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTAcatccatc                                                                                    >  1:404054/1‑114 (MQ=255)
                                                                                                          cTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATg                                          <  1:417462/150‑1 (MQ=255)
                                                                                                                   cTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCcaccac                                 >  2:1039562/1‑150 (MQ=255)
                                                                                                                    tGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCaccacc                                >  1:488041/1‑150 (MQ=255)
                                                                                                                     gagccATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCAccc                               >  1:2330452/4‑150 (MQ=255)
                                                                                                                     gATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTAt                                           <  2:332197/138‑1 (MQ=255)
                                                                                                                     gATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTAt                                           >  1:332197/1‑138 (MQ=255)
                                                                                                                      aTCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCAcccc                              <  1:19110/150‑1 (MQ=255)
                                                                                                                       tCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCg                             >  1:1661117/1‑150 (MQ=255)
                                                                                                                       tCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCg                             >  2:2023072/1‑150 (MQ=255)
                                                                                                                        ccATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTAcatccatc                                                                                    <  2:2081943/94‑1 (MQ=255)
                                                                                                                        ccATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTAcatccatc                                                                                    >  1:2081943/1‑94 (MQ=255)
                                                                                                                         cATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGAt                           >  2:1803227/1‑150 (MQ=255)
                                                                                                                          atgaCCCGCACGGCCCCCTTCGCGGGTTTACCCGCGAAGAGGTCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATg                          <  2:252511/147‑1 (MQ=255)
                                                                                                                          aTCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATcc                                              <  2:488011/130‑1 (MQ=255)
                                                                                                                          aTCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATcc                                              >  1:488011/1‑130 (MQ=255)
                                                                                                                           tCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGa                         >  2:1864172/1‑150 (MQ=255)
                                                                                                                            cACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGaa                        <  1:1108508/150‑1 (MQ=255)
                                                                                                                               cTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCg                     >  1:1950821/1‑150 (MQ=255)
                                                                                                                                 gCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTAcatcca                                                                                      <  1:841516/83‑1 (MQ=255)
                                                                                                                                 gCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTAcatcca                                                                                      >  2:841516/1‑83 (MQ=255)
                                                                                                                                       gCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCAccc                               >  1:838084/1‑132 (MQ=255)
                                                                                                                                       gCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCAccc                               <  2:838084/132‑1 (MQ=255)
                                                                                                                                       gCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGc             >  1:922125/1‑150 (MQ=255)
                                                                                                                                         ctctTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGGTAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCa           >  1:1996830/1‑150 (MQ=255)
                                                                                                                                            ttCGCGGGTAAACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATgc        >  2:562953/1‑150 (MQ=255)
                                                                                                                                                 gggTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATGCACCTg   >  1:107469/1‑150 (MQ=255)
                                                                                                                                                  ggTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATGCGCCTGt  <  1:1803227/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
TCAGCGCCCAGGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATGCGCCTGT  >  NC_002947/4529920‑4530215

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