Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,263,006 A→G 6.6% N165N (AAT→AAC lolB ← outer membrane lipoprotein LolB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,263,0060AG 6.6% 274.0 / 7.4 91N165N (AAT→AAClolBouter membrane lipoprotein LolB
Reads supporting (aligned to +/- strand):  ref base A (45/40);  new base G (3/3);  total (48/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

GCCACTGTGTCCGCATTATTTCACTATCCAGTTATCCATTTTTAACTTGATGCGTTGACCACCGTCGGTGAGTTCCATATTGGCTGGCATCGCAGGTTGCGTTTTGGTGTCATAACCACCATAAACAACCTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGGTCGTCCAGTTTGTAGTCGGTTGCATCACCCGGTAAACCTAAAATCCACTGGCGCAAGCTGTTGAGCGGAATTGGCATTCCGGTCAATTTGCCAATCATCTCTTC  >  NC_000913/1262862‑1263146
                                                                                                                                                |                                                                                                                                            
gCCACTGTGTCCGCATTATTTCACTATCCAGTTATCCATTTTTAACTTGATGCGTTGACCACCGTCGGTGAGTTCCATATTGGCTGGCATCGCAGGTTGCGTTTTGGTGTCATAACCACCATAAACAACCTTCCAGTTTTTGCCATTCTg                                                                                                                                         >  1:650362/1‑150 (MQ=255)
 ccACTGTGTCCGCATTATTTCACTATCCAGTTATCCATTTTTAACTTGATGCGTTGACCACCGTCGGTGAGTTCCATATTGGCTGGCATCGCAGGTTGCGTTTTGGTGTCATAACCACCATAAACAACCTTCCAGTTTTTGCCATTCTgg                                                                                                                                        >  1:1254758/1‑150 (MQ=255)
  cACTGTGTCCGCATTATTTCACTATCCAGTTATCCATTTTTAACTTGATGCGTTGACCACCGTCGGTGAGTTCCATATTGGCTGGCATCGCAGGTTGCGTTTTGGTGTCATAACCACCATAAACAACCTTCCAGTTTTTGCCATTCTGGc                                                                                                                                       <  1:1036623/150‑1 (MQ=255)
  cACTGTGTCCGCATTATTTCACTATCCAGTTATCCATTTTTAACTTGATGCGTTGACCACCGTCGGTGAGTTCCATATTGGCTGGCATCGCAGGTTGCGTTTTGGTGTCATAACCACCATAAACAACCTTCCAGTTTTTGCCATTCTGGc                                                                                                                                       >  2:1226981/1‑150 (MQ=255)
  cACTGTGTCCACATTATTTCACTATCCAGTTATCCATTTTTAACTTGATGCGTTGACCACCGTCGGTGAGTTCCATATTGGCTGGCATCGCAGGTTGCGTTTTGGTGTCATAACCACCATAAACAACCTTCCAGTTTTTGCCATTCTGGc                                                                                                                                       <  1:204319/150‑1 (MQ=255)
   aCTGTGTCCGCATTATTTCACTATCCAGTTATCCATTTTTAACTTGATGCGTTGACCACCGTCGGTGAGTTCCATATTGGCTGGCATCGCAGGTTGCGTTTTGGTGTCATAACCACCATAAACAACCTTCCAGTTTTTGCCATTCTGGCt                                                                                                                                      <  2:175205/150‑1 (MQ=255)
    cTGTGTCCGCATTATTTCACTATCCAGTTATCCATTTTTAACTTGATGCGTTGACCACCGTCGGTGAGTTCCATATTGGCTGGCATCGCAGGTTGCGTTTTGGTGTCATAACCACCATAAACAACCTTCCAGTTTTTGCCATTCTGGCTg                                                                                                                                     >  2:150650/1‑150 (MQ=255)
          ccGCATTATTTCACTATCCAGTTATCCATTTTTAACTTGATGCGTTGACCACCGTCGGTGAGTTCCATATTGGCTGGCATCGCAGGTTGCGTTTTGGTGTCATAACCACCATAAACAACCTTCCAGTTTTTGCCATTCTGGCTGTAGGTa                                                                                                                               >  2:1590769/1‑150 (MQ=255)
          ccGCATTATTTCACTATCCAGTTATCCATTTTTAACTTGATGCGTTGACCACCGTCGGTGAGTTCCATATTGGCTGGCATCGCAGGTTGCGTTTTGGTGTCATAACCACCATAAACAACCTTCCAGTTTTTGCCATTCTGGCTGTAGGTa                                                                                                                               >  2:1608122/1‑150 (MQ=255)
                tattTCACTATCCAGTTATCCATTTTTAACTTGATGCGTTGACCACCGTCGGTGAGTTCCATATTGGCTGGCATCGCAGGTTGCGTTTTGGTGTCATAACCACCATAAACAACCTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCg                                                                                                                         >  2:878265/1‑150 (MQ=255)
                  tttCACTATCCAGTTATCCATTTTTAACTTGATGCGTTGACCACCGTCGGTGAGTTCCATATTGGCTGGCATCGCAGGTTGCGTTTTGGTGTCATAACCACCATAAACAACCTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCt                                                                                                                       >  1:1283200/1‑150 (MQ=255)
                  tttCACTATCCAGTTATCCATTTTTAACTTGATGCGTTGACCACCGTCGGTGAGTTCCATATTGGCTGGCATCGCAGGTTGCGTTTTGGTGTCATAACCACCATAAACAACCTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCt                                                                                                                       >  1:1277137/1‑150 (MQ=255)
                   ttCACTATCCAGTTATCCATTTTTAACTTGATGCGTTGACCACCGTCGGTGAGTTCCATATTGGCTGGCATCGCAGGTTGCGTTTTGGTGTCATAACCACCATAAACAACCTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTc                                                                                                                      <  1:150650/150‑1 (MQ=255)
                    tCACTATCCAGTTATCCATTTTTAACTTGATGCGTTGACCACCGTCGGTGAGTTCCATATTGGCTGGCATCGCAGGTTGCGTTTTGGTGTCATAACCACCATAAACAACCTTCCAGTTTTTGCCGTTCTGGCTGTAGGTAATTTCGCTCa                                                                                                                     >  1:845040/1‑150 (MQ=255)
                    tCACTATCCAGTTATCCATTTTTAACTTGATGCGTTGACCACCGTCGGTGAGTTCCATATTGGCTGGCATCGCAGGTTGCGTTTTGGTGTCATAACCACCATAAACAACCTTCCAGTTTTTGCCGTTCTGGCTGTAGGTAATTTCGCTCa                                                                                                                     >  1:309200/1‑150 (MQ=255)
                    tCACTATCCAGTTATCCATTTTTAACTTGATGCGTTGACCACCGTCGGTGAGTTCCATATTGGCTGGCATCGCAGGTTGCGTTTTGGTGTCATAACCACCATAAACAACCTTCCAGTTTTTGCCGTTCTGGCTGTAGGTAATTTCGCTCa                                                                                                                     >  1:850664/1‑150 (MQ=255)
                       cTATCCAGTTATCCATTTTTAACTTGATGCGTTGACCACCGTCGGTGAGTTCCATATTGGCTGGCATCGCAGGTTGCGTTTTGGTGTCATAACCACCATAAACAACCTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGc                                                                                                                  <  1:1302400/150‑1 (MQ=255)
                       cTATCCAGTTATCCATTTTTAACTTGATGCGTTGACCACCGTCGGTGAGTTCCATATTGGCTGGCATCGCAGGTTGCGTTTTGGTGTCATAACCACCATAAACAACCTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGc                                                                                                                  >  1:325985/1‑150 (MQ=255)
                         aTCCAGTTATCCATTTTTAACTTGATGCGTTGACCACCGTCGGTGAGTTCCATATTGGCTGGCATCGCAGGTTGCGTTTTGGTGTCATAACCACCATAAACAACCTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCt                                                                                                                       >  1:1602534/1‑143 (MQ=255)
                         aTCCAGTTATCCATTTTTAACTTGATGCGTTGACCACCGTCGGTGAGTTCCATATTGGCTGGCATCGCAGGTTGCGTTTTGGTGTCATAACCACCATAAACAACCTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCt                                                                                                                       <  2:1602534/143‑1 (MQ=255)
                            cAGTTATCCATTTTTAACTTGATGCGTTGACCACCGTCGGTGAGTTCCATATTGGCTGGCATCGCAGGTTGCGTTTTGGTGTCATAACCACCATAATCAACCTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTAc                                                                                                             >  1:51473/1‑150 (MQ=255)
                             aGTTATCCATTTTTAGCTTGATGCGTTGACCACCGTCGGTGAGTTCCATATTGGCTGGCATCGCAGGTTGCGTTTTGGTGTCATAACCACCATAAACAACCTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACt                                                                                                            <  1:1307955/150‑1 (MQ=255)
                             aGTTATCCATTTTTAACTTGATGCGTTGACCACCGTCGGTGAGTTCCATATTGGCTGGCATCGCAGGTTGCGTTTTGGTGTCATAACCACCATAAACAACCTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACt                                                                                                            >  1:1621503/1‑150 (MQ=255)
                                tATCCATTTTTAACTTGATGCGTTGACCACCGTCGGTGAGTTCCATATTGGCTGGCATCGCAGGTTGCGTTTTGGTGTCATAACCACCATAATCAACCTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGgt                                                                                                         <  2:51473/150‑1 (MQ=255)
                                tATCCATTTTTAACTTGATGCGTTGACCACCGTCGGTGAGTTCCATATTGGCTGGCATCGCAGGTTGCGTTTTGGTGTCATAACCACCATAAACAACCTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGgt                                                                                                         <  1:998509/150‑1 (MQ=255)
                                 aTCCATTTTTAACTTGATGCGTTGACCACCGTCGGTGAGTTCCATATTGGCTGGCATCGCAGGTTGCGTTTTGGTGTCATAACCACCATAAACAACCTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGgtc                                                                                                        >  2:92759/1‑150 (MQ=255)
                                   ccATTTTTAACTTGATGCGTTGACCACCGTCGGTGAGTTCCATATTGGCTGGCATCGCAGGTTGCGTTTTGGTGTCATAACCACCATAAACAACCTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGgtcgt                                                                                                      >  2:576201/1‑150 (MQ=255)
                                   ccATTTTTAACTTGATGCGTTGACCACCGTCGGTGAGTTCCATATTGGCTGGCATCGCAGGTTGCGTTTTGGTGTCATAACCACCATAAACAACCTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGgtcgt                                                                                                      >  1:1127367/1‑150 (MQ=255)
                                    cATTTTTAACTTGATGCGTTGACCACCGTCGGTGAGTTCCATATTGGCTGGCATCGCAGGTTGCGTTTTGGTGTCATAACCACCATAAACAACCTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGgtcgtc                                                                                                     <  2:139763/150‑1 (MQ=255)
                                          tAACTTGATGCGTTGACCACCGTCGGTGAGTTCCATATTGGCTGGCATCGCAGGTTGCGTTTTGGTGTCATAACCACCATAAACAACCTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGGTCGTCCAGttt                                                                                               >  1:751488/1‑150 (MQ=255)
                                            aCTTGATGCGTTGACCACCGTCGGTGAGTTCCATATTGGCTGGCATCGCAGGTTGCGTTTTGGTGTCATAACCACCATAAACAACCTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGGTCGTCCAGTTTGt                                                                                             <  2:650362/150‑1 (MQ=255)
                                            aCTTGATGCGTTGACCACCGTCGGTGAGTTCCATATTGGCTGGCATCGCAGGTTGCGTTTTGGTGTCATAACCACCATAAACAACCTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGGTCGTCCAGTTTGt                                                                                             <  1:878265/150‑1 (MQ=255)
                                            aCTTGATGCGTTGACCACCGTCGGTGAGTTCCATATTGGCTGGCATCGCAGGTTGCGTTTTGGTGTCATAACCACCATAAACAACCTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGGTCGTCCAGTTTGt                                                                                             <  2:187689/150‑1 (MQ=255)
                                            aCTTGATGCGTTGACCACCGTCGGTGAGTTCCATATTGGCTGGCATCGCAGGTTGCGTTTTGGTGGCATAACCACCATAAACAACCTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGGGCGTCCAGTTTGt                                                                                             <  2:1153618/150‑1 (MQ=255)
                                             cTTGATGCGTTGACCACCGTCGGTGAGTTCCATATTGGCTGGCATCGCAGGTTGCGTTTTGGTGTCATAACCACCATAAACAACCTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGGTCGTCCAGTTTGTa                                                                                            >  2:251464/1‑150 (MQ=255)
                                              ttGATGCGTTGACCACCGTCGGTGAGTTCCATATTGGCTGGCATCGCAGGTTGCGTTTTGGTGTCATAACCACCATAAACAACCTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGGTCGTCCAGTTTGTAg                                                                                           <  1:1449094/150‑1 (MQ=255)
                                                  tGCGTTGCCCACCGTCGGTGAGTTCCATATTGGCTGGCATCGCAGGTTGCGTTTTGGTGTCATAACCCCCATAAACACCCTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGGTCGTCCAGTTTGTAGTCgg                                                                                       >  1:1370685/1‑150 (MQ=255)
                                                    cGTTGACCACCGTCGGTGAGTTCCATATTGGCTGGCATCGCAGGTTGCGTTTTGGTGTCATAACCACCATAAACAACCTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGGTCGTCCAGTTTGTAGTCGGtt                                                                                     >  1:646977/1‑150 (MQ=255)
                                                     gTTGACCACCGTCGGTGAGTTCCATATTGGCTGGCATCGCAGGTTGCGTTTTGGTGTCATAACCACCATAAACAACCTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGGTCGTCCAGTTTGTAGTCGGTTg                                                                                    <  2:656744/150‑1 (MQ=255)
                                                      ttGACCACCGTCGGTGAGTTCCATATTGGCTGGCATCGCAGGTTGCGTTTTGGTGTCATAACCACCATAAACAACCTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGGTCGTCCAGTTTGTAGTCGGTTGc                                                                                   <  1:370282/150‑1 (MQ=255)
                                                        gACCACCGTCGGTGAGTTCCATATTGGCTGGCATCGCAGGTTGCGTTTTGGTGTCATAACCACCATAAACAACCTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGGTCGTCCAGTTTGTAGTCGGTTGCAt                                                                                 <  1:514822/150‑1 (MQ=255)
                                                          ccaccGTCGGTGAGTTCCATATTGGCTGGCATCGCAGGTTGCGTTTTGGTGTCATAACCACCATAAACAACCTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGGTCGTCCAGTTTGTAGTCGGTTGCATCa                                                                               >  1:654459/1‑150 (MQ=255)
                                                          ccaccGTCGGTGAGTTCCATATTGGCTGGCATCGCAGGTTGCGTTTTGGTGTCATAACCACCATAAACAACCTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGGTCGTCCAGTTTGTAGTCGGTTGCATCa                                                                               <  2:341136/150‑1 (MQ=255)
                                                             ccGTCGGTGAGTTCCATATTTGCTGGCATCGCAGGGTGCGTTTTGTGGTGTTAACCACCATAAACAACCTTCCAGTTTTTGCCATTATGGCTGTAGTTAATTTCGCTCAGGCGGTACTGGTCGTCCTGTTTGTAGTCGGTTGCATCAccc                                                                            <  2:779402/150‑1 (MQ=255)
                                                             ccGTCGGTGAGTTCCATATTGGCTGGCATCGCAGGTTGCGTTTTGGTGTCATAACCACCATAAACAACCTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGGTCGTCCAGTTTGTAGTCGGTTGCATCAccc                                                                            <  2:1521051/150‑1 (MQ=255)
                                                              cGTCGGTGAGTTCCATATTGGCTGGCATCGCAGGTTGCGTTTTGGTGTCATAACCACCATAAACAACCTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGGTCGTCCAGTTTGTAGTCGGTTGCATCACCCg                                                                           <  1:1534538/150‑1 (MQ=255)
                                                                 cGGTGAGTTCCATATTGGCTGGCATCGCAGGTTGCGTTTTGGTGTCATAACCACCATAAACAACCTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGGTCGTCCAGTTTGTAGTCGGTTGCATCACCCGGTa                                                                        <  2:982725/150‑1 (MQ=255)
                                                                   gTGAGTTCCATATTGGCTGGCATCGCAGGTTGCGTTTTGGTGTCATAACCACCATAAACAACCTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGGTCGTCCAGTTTGTAGTCGGTTGCATCACCCGGTaaa                                                                      >  2:1026528/1‑150 (MQ=255)
                                                                     gAGTTCCATATTGGCTGGCATCGCAGGTTGCGTTTTGGTGTCATAACCACCATAAACAACCTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGGTCGTCCAGTTTGTAGTCGGTTGCATCACCCGGTAAAcc                                                                    <  1:908152/150‑1 (MQ=255)
                                                                     gAGTTCCATATTGGCTGGCATCGCAGGTTGCGTTTTGGTGTCATAACCACCATAAACAACCTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGGTCGTCCAGTTTGTAGTCGGTTGCATCACCCGGTAAAcc                                                                    <  1:576201/150‑1 (MQ=255)
                                                                      aGTTCCATATTGGCTGGCATCGCAGGTTGCGTTTTGGTGTCATAACCACCATAAACAACCTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGGTCGTCCAGTTTGTAGTCGGTTGCATCACCCGGTAAACCt                                                                   >  2:484598/1‑150 (MQ=255)
                                                                         tCCATATTGGCTGGCATCGCAGGTTGCGTTTTGGTGTCATAACCACCATAAACAACCTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGGTCGTCCAGTTTGTAGTCGGTTGCATCACCCGGTAAACCTaaa                                                                >  1:462367/1‑150 (MQ=255)
                                                                          ccATATTGGCTGGCATCGCAGGTTGCGTTTTGGTGTCATAACCACCATAAACAACCTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGGTCGTCCAGTTTGTAGTCGGTGGCATCACCCGGTAAACCTaaaa                                                               <  2:655401/150‑1 (MQ=255)
                                                                            atatTGGCTGGCATCGCAGGTTGCGTTTTGGTGTCATAACCACCATAAACAACCTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGGTCGTCCAGTTTGTAGTCGGTTGCATCACCCGGTAAACCTAAAATc                                                             <  2:1596588/150‑1 (MQ=255)
                                                                              atTGGCTGGCATCGCAGGTTGCGTTTTGGTGTCATAACCACCATAAACAACCTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGGTCGTCCAGTTTGTAGTCGGTTGCATCACCCGGTAAACCTAAAATCCa                                                           >  2:850115/1‑150 (MQ=255)
                                                                                tggctggcATCGCAGGTTGCGTTTTGGTGTCATAACCACCATAAACAACCTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGGTCGTCCAGTTTGTAGTCGGTTGCATCACCCGGTAAACCTAAAATCCACt                                                         <  1:313542/150‑1 (MQ=255)
                                                                                  gctggcATCGCAGGTTGCGTTTTGGTGTCATAACCACCATAAACAACCTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGGTCGTCCAGTTTGTAGTCGGTTGCATCACCCGGTAAACCTAAAATCCACTgg                                                       >  2:201665/1‑150 (MQ=255)
                                                                                    tggcATCGCAGGTTGCGTTTTGGTGTCATAACCACCATAAACAACCTTCCAGTTTTTGCCGTTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGGTCGTCCAGTTTGTAGTCGGTTGCATCACCCGGTAAACCTAAAATCCACTGgcg                                                     <  2:845040/150‑1 (MQ=255)
                                                                                    tggcATCGCAGGTTGCGTTTTGGTGTCATAACCACCATAAACAACCTTCCAGTTTTTGCCGTTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGGTCGTCCAGTTTGTAGTCGGTTGCATCACCCGGTAAACCTAAAATCCACTGgcg                                                     <  2:850664/150‑1 (MQ=255)
                                                                                    tggcATCGCAGGTTGCGTTTTGGTGTCATAACCACCATAAACAACCTTCCAGTTTTTGCCGTTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGGTCGTCCAGTTTGTAGTCGGTTGCATCACCCGGTAAACCTAAAATCCACTGgcg                                                     <  2:309200/150‑1 (MQ=255)
                                                                                         tCGCAGGTTGCGTTTTGGTGTCATAACCACCATAAACAACCTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGGTCGTCCAGTTTGTAGTCGGTTGCATCACCCGGTAAACCTAAAATCCACTGGCGCAAGc                                                >  2:1579048/1‑150 (MQ=255)
                                                                                         tCGCAGGTTGCGTTTTGGTGTCATAACCACCATAAACAACCTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGGTCGTCCAGTTTGTAGTCGGTTGCATCACCCGGTAAACCTAAAATCCACTGGCGCAAGc                                                >  2:823245/1‑150 (MQ=255)
                                                                                            cAGGTTGCGTTTTGGTGTCATAACCACCATAAACAACCTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGGTCGTCCAGTTTGTAGTCGGTTGCATCACCCGGTAAACCTAAAATCCACTGGCGCAAGCTGt                                             <  2:1127367/150‑1 (MQ=255)
                                                                                            cAGGTTGCGTTTTGGTGTCATAACCACCATAAACAACCTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGGTCGTCCAGTTTGTAGTCGGTTGCATCACCCGGTAAACCTAAAATCCACTGGCGCAAGCTGt                                             <  1:1608122/150‑1 (MQ=255)
                                                                                              ggTTGCGTTTTGGTGTCATAACCACCATAAACAACCTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGGTCGTCCAGTTTGTAGTCGGTTGCATCACCCGGTAAACCTAAAATCCACTGGCGCAAGCTGTTg                                           <  1:1579048/150‑1 (MQ=255)
                                                                                               gTTGCGTTTTGGTGTCATAACCACCATAAACAACCTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGGTCGTCCAGTTTGTAGTCGGTTGCATCACCCGGTAAACCTAAAATCCACTGGCGCAAGCTGTTGa                                          >  1:784639/1‑150 (MQ=255)
                                                                                                   cGTTTTGGTGTCATAACCACCATAAACAACCTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGGTCGTCCAGTTTGTAGTCGGTTGCATCACCCGGTAAACCTAAAATCCACTGGCGCAAGCTGTTGAGCgg                                      >  2:806895/1‑150 (MQ=255)
                                                                                                        tGGTGTCATAACCACCATAAACAACCTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGGTCGTCCAGTTTGTAGTCGGTTGCATCACCCGGTAAACCTAAAATCCACTGGCGCAAGCTGTTGAGCGGAATTg                                 <  2:646977/150‑1 (MQ=255)
                                                                                                         ggTGTCATAACCACCATAAACAACCTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGGTCGTCCAGTTTGTAGTCGGTTGCATCACCCGGTAAACCTAAAATCCACTGGCGCAAGCTGTTGAGCGGAATTgg                                <  1:1226981/150‑1 (MQ=255)
                                                                                                            gtCATAACCACCATAAACAACCTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGGTCGTCCAGTTTGTAGTCGGTTGCATCACCCGGTAAACCTAAAATCCACTGGCGCAAGCTGTTGAGCGGAATTGGCAt                             >  2:484874/1‑150 (MQ=255)
                                                                                                            gtCATAACCACCATAAACAACCTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGGTCGTCCAGTTTGTAGTCGGTTGCATCACCCGGTAAACCTAAAATCCACTGGCGCAAGCTGTTGAGCGGAATTGGCAt                             >  2:402085/1‑150 (MQ=255)
                                                                                                            gtCATAACCACCATAAACAACCTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGGTCGTCCAGTTTGTAGTCGGTTGCATCACCCGGTAAACCTAAAATCCACTGGCGCAAGCTGTTGAGCGGAATTGGCAt                             >  2:261133/1‑150 (MQ=255)
                                                                                                               aTAACCACCATAAACAACCTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGGTCGTCCAGTTTGTAGTCGGTTGCATCACCCGGTAAACCTAAAATCCACTGGCGCAAGCTGTTGAGCGGAATTGGCATTcc                          >  2:642690/1‑150 (MQ=255)
                                                                                                                tAACCACCATAAACAACCTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGGTCGTCCAGTTTGTAGTCGGTTGCATCACCCGGTAAACCTAAAATCCACTGGCGCAAGCTGTTGAGCGGAATTGGCATTCCg                         >  1:1150680/1‑150 (MQ=255)
                                                                                                                 aaCCACCATAAACAACCTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGGTCGTCCAGTTTGTAGTCGGTTGCATCACCCGGTAAACCTAAAATCCACTGGCGCAAGCTGTTGAGCGGAATTGGCATTCCgg                        <  2:1621094/150‑1 (MQ=255)
                                                                                                                   ccaccaTAAACAACCTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGGTCGTCCAGTTTGTAGTCGGTTGCATCACCCGGTAAACCTAAAATCCACTGGCGCAAGCTGTTGAGCGGAATTGGCATTCCGGTc                      <  2:1597468/150‑1 (MQ=255)
                                                                                                                      ccaTAAACAACCTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGGTCGTCCAGTTTGTAGTCGGTTGCATCACCCGGTAAACCTAAAATCCACTGGCGCAAGCTGTTGAGCGGAATTGGCATTCCGGTCAAt                   >  1:895706/1‑150 (MQ=255)
                                                                                                                               acCTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGGTCGTCCAGTTTGTAGTCGGTTGCATCACCCGGTAAACCTAAAATCCACTGGCGCAAGCTGTTGAGCGGAATTGGCATTCCGGTCAATTTGCCAatc          >  1:1526968/1‑150 (MQ=255)
                                                                                                                                ccTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGGTCGTCCAGTTTGTAGTCGGTTGCATCACCCGGTAAACCTAAAATCCACTGGCGCAAGc                                                >  2:611001/1‑111 (MQ=255)
                                                                                                                                ccTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGGTCGTCCAGTTTGTAGTCGGTTGCATCACCCGGTAAACCTAAAATCCACTGGCGCAAGc                                                <  1:611001/111‑1 (MQ=255)
                                                                                                                                ccTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGGTCGTCCAGTTTGTAGTCGGTTGCATCACCCGGTAAACCTAAAATCCACTGGCGCAAGCTGTTGAGCGGAATTGGCATTCCGGTCAATTTGCCAatca         <  2:751488/150‑1 (MQ=255)
                                                                                                                                ccTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGGTCGTCCAGTTTGTAGTCGGTTGCATCACCCGGTAAACCTAAAATCCACTGGCGCAAGCTGTTGAGCGGAATTGGCATTCCGGTCAATTTGCCAatca         >  1:358584/1‑150 (MQ=255)
                                                                                                                                 cTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGGTCGTCCAGTTTGTAGTCGGTTGCATCACCCGGTAAACCTAAAATCCACTGGCGCAAGCTGTTGAGCGGAATTGGCATTCCGGTCAATTTGCCAatcat        >  1:1039481/1‑150 (MQ=255)
                                                                                                                                  ttCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGGTCGTCCAGTTTGTAGTCGGTTGCATCACCCGGTAAACCTAAAATCCACTGGCGCAAGCTGTTGAGCGGAATTGGCATTCCGGTCAATTTGCCAatcatc       >  1:1184285/1‑150 (MQ=255)
                                                                                                                                  ttCCAGTTTTTGCCATTCTGGCTGTAGCTAATTTCGCTCAGGCCGTACTCGTCGTCCAGTTTGTAGTCGGTTGCATCACCCGCTAAACCTAAAATCCACTGGCGCAATCTGTTGAGCGGAATTGGCATTCCGGTCAATTTGCCAatcatc       >  1:117757/1‑150 (MQ=255)
                                                                                                                                   tCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGGTCGTCCAGTTTGTAGTCGGTTGCATCACCCGGTAAACCTAAAATCCACTGGCGCAAGCTGTTGAGCGGAATTg                                 <  1:552027/123‑1 (MQ=255)
                                                                                                                                   tCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGGTCGTCCAGTTTGTAGTCGGTTGCATCACCCGGTAAACCTAAAATCCACTGGCGCAAGCTGTTGAGCGGAATTg                                 >  2:552027/1‑123 (MQ=255)
                                                                                                                                   tCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGGTCGTCCAGTTTGTAGTCGGTTGCATCACCCGGTAAACCTAAAATCCACTGGCGCAAGCTGTTGAGCGGAATTGGCATTCCGGTCAATTTGCCAATCAtct      <  2:325985/150‑1 (MQ=255)
                                                                                                                                   tCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGGTCGTCCAGTTTGTAGTCGGTTGCATCACCCGGTAAACCTAAAATCCACTGGCGCAAGCTGTTGAGCGGAATTGGCATTCCGGTCAATTTGCCAATCAtct      <  1:642690/150‑1 (MQ=255)
                                                                                                                                    ccAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGGTCGTCCAGTTTGTAGTCGGTTGCATCACCCGGTAAACCTAAAATCCACTGGCGCAAGCTGTTGAGCGGAATTGGCATTCCGGTCAATTTGCCAATCAtctc     >  2:1311565/1‑150 (MQ=255)
                                                                                                                                       gTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGGTCGTCCAGTTTGTAGTCGGTTGCATCACCCGGTAAACCTAAAATCCACTGGCGCAAGCTGTTGAGCGGAATTGGCATTCCGGTCAATTTGCCAATCATCTCTTc  >  1:344826/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                            
GCCACTGTGTCCGCATTATTTCACTATCCAGTTATCCATTTTTAACTTGATGCGTTGACCACCGTCGGTGAGTTCCATATTGGCTGGCATCGCAGGTTGCGTTTTGGTGTCATAACCACCATAAACAACCTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGGTCGTCCAGTTTGTAGTCGGTTGCATCACCCGGTAAACCTAAAATCCACTGGCGCAAGCTGTTGAGCGGAATTGGCATTCCGGTCAATTTGCCAATCATCTCTTC  >  NC_000913/1262862‑1263146

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