Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,368,528 C→A R150R (CGC→CGA pspA → phage shock protein A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,368,5280CA98.8% 290.8 / NA 83R150R (CGC→CGApspAphage shock protein A
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (37/45);  total (38/45)
Fisher's exact test for biased strand distribution p-value = 4.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATGAAGTGACGCTGGTGGACGATACGCTGGCACGCATGAAGAAAGAGATTGGCGAGCTGGAAAACAAATTGAGCGAAACACGCGCTCGCCAGCAGGCATTGATGTTACGTCATCAGGCGGCAAACTCGTCGCGCGATGTGCGTCGTCAGCTGGACAGTGGCAAACTGGATGAAGCAATGGCTCGTTTCGAATCTTTCGAACGTCGTATTGACCAGATGGAAGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGAT  >  NC_000913/1368394‑1368663
                                                                                                                                      |                                                                                                                                       
cATGAAGTGACGCTGGTGGACGATACGCTGGCACGCATGAAGAAAGAGATTGGCGAGCTGGAAAACAAATTGAGCGAAACACGCGCTCGCCAGCAGGCATTGATGTTACGTCATCAGGCGGCAAACTCGTCGCGAGAtgt                                                                                                                                    <  1:1210673/140‑1 (MQ=255)
  tGAAGTGACGCTGGTGGACGATACGCTGGCACGCATGAAGAAAGAGATTGGCGAGCTGGAAAACAAATTGAGCGAAACACGCGCTCGCCAGCAGGCATTGATGTTACGTCATCAGGCGGCAAACTCGTCGCGAGATGTGc                                                                                                                                  >  2:792908/1‑140 (MQ=255)
     aGTGACGCTGGTGGACGATACGCTGGCACGCATGAAGAAAGAGATTGGCGAGCTGGAAAACAAATTGAGCGAAACACGCGCTCGCCAGCAGGCATTGATGTTACGTCATCAGGCGGCAAACTCGTCGCGAGATGTGcgtc                                                                                                                               <  2:1679943/140‑1 (MQ=255)
     aGTGACGCTGGTGGACGATACGCTGGCACGCATGAAGAAAGAGATTGGCGAGCTGGAAAACAAATTGAGCGAAACACGCGCTCGCCAGCAGGCATTGATGTTACGTCATCAGGCGGCAAACTCGTCGCGAGATGTGcgtc                                                                                                                               <  1:792908/140‑1 (MQ=255)
     aGTGACGCTGGTGGACGATACGCTGGCACGCATGAAGAAAGAGATTGGCGAGCTGGAAAACAAATTGAGCGAAACACGCGCTCGCCAGCAGGCATTGATGTTACGTCATCAGGCGGCAAACTCGTCGCGAGATGTGcgtc                                                                                                                               <  1:232329/140‑1 (MQ=255)
     aGTGACGCTGGTGGACGATACGCTGGCACGCATGAAGAAAGAGATTGGCGAGCTGGAAAACAAATTGAGCGAAACACGCGCTCGCCAGCAGGCATTGATGTTACGTCATCAGGCGGCAAACTCGTCGCGAGATGTGcgtc                                                                                                                               <  1:1419828/140‑1 (MQ=255)
       tGACGCTGGTGGACGATACGCTGGCACGCATGAAGAAAGAGATTGGCGAGCTGGAAAACAAATTGAGCGAAACACGCGCTCGCCAGCAGGCATTGATGTTACGTCATCAGGCGGCAAACTCGTCGCGAGATGTGcgtcgt                                                                                                                             >  2:380616/1‑140 (MQ=255)
       tGACGCTGGTGGACGATACGCTGGCACGCATGAAGAAAGAGATTGGCGAGCTGGAAAACAAATTGAGCGAAACACGCGCTCGCCAGCAGGCATTGATGTTACGTCATCAGGCGGCAAACTCGTCGCGAGATGTGcgtcgt                                                                                                                             >  2:1002260/1‑140 (MQ=255)
               gtggACGATACGCTGGCACGCATGAAGAAAGAGATTGGCGAGCTGGAAAACAAATTGAGCGAAACACGCGCTCGCCAGCAGGCATTGATGTTACGTCATCAGGCGGCAAACTCGTCGCGAGATGTGCGTCGTCAGCTGGa                                                                                                                     >  2:988648/1‑140 (MQ=255)
               gtggACGATACGCTGGCACGCATGAAGAAAGAGATTGGCGAGCTGGAAAACAAATTGAGCGAAACACGCGCTCGCCAGCAGGCATTGATGTTACGTCATCAGGCGGCAAACTCGTCGCGAGATGTGCGTCGTCAGCTGGa                                                                                                                     >  1:1323469/1‑140 (MQ=255)
                tggACGATACGCTGGCACGCATGAAGAAAGAGATTGGCGAGCTGGAAAACAAATTGAGCGAAACACGCGCTCGCCAGCAGGCATTGATGTTACGTCATCAGGCGGCAAACTCGTCGCGAGATGTGCGTCGTCAGCTGGAc                                                                                                                    <  1:1002260/140‑1 (MQ=255)
                 ggACGATACGCTGGCACGCATGAAGAAAGAGATTGGCGAGCTGGAAAACAAATTGAGCGAAACACGCGCTCGCCAGCAGGCATTGATGTTACGTCATCAGGCGGCAAACTCGTCGCGAGATGTGCGTCGTCAGCTGGACa                                                                                                                   <  2:1287029/140‑1 (MQ=255)
                    cGATACGCTGGCACGCATGAAGAAAGAGATTGGCGAGCTGGAAAACAAATTGAGCGAAACACGCGCTCGCCAGCAGGCATTGATGTTACGTCATCAGGCGGCAAACTCGTCGCGAGATGTGCGTCGTCAGCTGGACAGTg                                                                                                                >  2:1635350/1‑140 (MQ=255)
                     gATACGCTGGCACGCATGAAGAAAGAGATTGGCGAGCTGGAAAACAAATTGAGCGAAACACGCGCTCGCCAGCAGGCATTGATGTTACGTCATCAGGCGGCAAACTCGTCGCGAGATGTGCGTCGTCAGCTGGACAGTgg                                                                                                               <  1:881804/140‑1 (MQ=255)
                       tACGCTGGCACGCATGAAGAAAGAGATTGGCGAGCTGGAAAACAAATTGAGCGAAACACGCGCTCGCCAGCAGGCATTGATGTTACGTCATCAGGCGGCAAACTCGTCGCGAGATGTGCGTCGTCAGCTGGACAGTGGCa                                                                                                             >  1:1185003/1‑140 (MQ=255)
                       tACGCTGGCACGCATGAAGAAAGAGATTGGCGAGCTGGAAAACAAATTGAGCGAAACACGCGCTCGCCAGCAGGCATTGATGTTACGTCATCAGGCGGCAAACTCGTCGCGAGATGTGCGTCGTCAGCTGGACAGTGGCa                                                                                                             >  1:1235656/1‑140 (MQ=255)
                       tACGCTGGCACGCATGAAGAAAGAGATTGGCGAGCTGGAAAACAAATTGAGCGAAACACGCGCTCGCCAGCAGGCATTGATGTTACGTCATCAGGCGGCAAACTCGTCGCGAGATGTGCGTCGTCAGCTGGACAGTGGCa                                                                                                             <  2:1405783/140‑1 (MQ=255)
                       tACGCTGGCACGCATGAAGAAAGAGATTGGCGAGCTGGAAAACAAATTGAGCGAAACACGCGCTCGCCAGCAGGCATTGATGTTACGTCATCAGGCGGCAAACTCGTCGCGAGATGTGCGTCGTCAGCTGGACAGTGGCa                                                                                                             >  1:1620459/1‑140 (MQ=255)
                         cGCTGGCACGCATGAAGAAAGAGATTGGCGAGCTGGAAAACAAATTGAGCGAAACACGCGCTCGCCAGCAGGCATTGATGTTACGTCATCAGGCGGCAAACTCGTCGCGAGATGTGCGTCGTCAGCTGGACAGTGGCaaa                                                                                                           >  2:404976/1‑140 (MQ=255)
                         cGCTGGCACGCATGAAGAAAGAGATTGGCGAGCTGGAAAACAAATTGAGCGAAACACGCGCTCGCCAGCAGGCATTGATGTTACGTCATCAGGCGGCAAACTCGTCGCGAGATGTGCGTCGTCAGCTGGACAGTGGCaaa                                                                                                           >  2:1507498/1‑140 (MQ=255)
                         cGCTGGCACGCATGAAGAAAGAGATTGGCGAGCTGGAAAACAAATTGAGCGAAACACGCGCTCGCCAGCAGGCATTGATGTTACGTCATCAGGCGGCAAACTCGTCGCGAGATGTGCGTCGTCAGCTGGACAGTGGCaaa                                                                                                           >  1:1434323/1‑140 (MQ=255)
                         cGCTGGCACGCATGAAGAAAGAGATTGGCGAGCTGGAAAACAAATTGAGCGAAACACGCGCTCGCCAGCAGGCATTGATGTTACGTCATCAGGCGGCAAACTCGTCGCGAGATGTGCGTCGTCAGCTGGACAGTGGCaaa                                                                                                           <  1:1438354/140‑1 (MQ=255)
                         cGCTGGCACGCATGAAGAAAGAGATTGGCGAGCTGGAAAACAAATTGAGCGAAACACGCGCTCGCCAGCAGGCATTGATGTTACGTCATCAGGCGGCAAACTCGTCGCGAGATGTGCGTCGTCAGCTGGACAGTGGCaaa                                                                                                           <  1:179980/140‑1 (MQ=255)
                         cGCTGGCACGCATGAAGAAAGAGATTGGCGAGCTGGAAAACAAATTGAGCGAAACACGCGCTCGCCAGCAGGCATTGATGTTACGTCATCAGGCGGCAAACTCGTCACGAGATGTGCGTCGTCAGCTGGACAGTGGCaaa                                                                                                           <  2:1235656/140‑1 (MQ=255)
                              gCACGCATGAAGAAAGAGATTGGCGAGCTGGAAAACAAATTGAGCGAAACACGCGCTCGCCAGCAGGCATTGATGTTACGTCATCAGGCGGCAAACTCGTCGCGAGATGTGCGTCGTCAGCTGGACAGTGGCAAACTGGa                                                                                                      <  1:1313110/140‑1 (MQ=255)
                              gCACGCATGAAGAAAGAGATTGGCGAGCTGGAAAACAAATTGAGCGAAACACGCGCTCGCCAGCAGGCATTGATGTTACGTCATCAGGCGGCAAACTCGTCGCGAGATGTGCGTCGTCAGCTGGACAGTGGCAAACTGGa                                                                                                      <  2:1586726/140‑1 (MQ=255)
                                 cGCATGAAGAAAGAGATTGGCGAGCTGGAAAACAAATTGAGCGAAACACGCGCTCGCCAGCAGGCATTGATGTTACGTCATCAGGCGGCAAACTCGTCGCGAGATGTGCGTCGTCAGCTGGACAGTGGCAAACTGGATGa                                                                                                   >  1:768470/1‑140 (MQ=255)
                                 cGCATGAAGAAAGAGATTGGCGAGCTGGAAAACAAATTGAGCGAAACACGCGCTCGCCAGCAGGCATTGATGTTACGTCATCAGGCGGCAAACTCGTCGCGAGATGTGCGTCGTCAGCTGGACAGTGGCAAACTGGATGa                                                                                                   >  2:1129355/1‑140 (MQ=255)
                                 cGCATGAAGAAAGAGATTGGCGAGCTGGAAAACAAATTGAGCGAAACACGCGCTCGCCAGCAGGCATTGATGTTACGTCATCAGGCGGCAAACTCGTCGCGAGATGTGCGTCGTCAGCTGGACAGTGGCAAACTGGATGa                                                                                                   <  2:1323469/140‑1 (MQ=255)
                                 cGCATGAAGAAAGAGATTGGCGAGCTGGAAAACAAATTGAGCGAAACACGCGCTCGCCAGCAGGCATTGATGTTACGTCATCAGGCGGCAAACTCGTCGCGAGATGTGCGTCGTCAGCTGGACAGTGGCAAACTGGATGa                                                                                                   <  1:1507498/140‑1 (MQ=255)
                                   cATGAAGAAAGAGATTGGCGAGCTGGAAAACAAATTGAGCGAAACACGCGCTCGCCAGCAGGCATTGATGTTACGTCATCAGGCGGCAAACTCGTCGCGAGATGTGCGTCGTCAGCTGGACAGTGGCaaa                                                                                                           <  2:1024033/130‑1 (MQ=255)
                                   cATGAAGAAAGAGATTGGCGAGCTGGAAAACAAATTGAGCGAAACACGCGCTCGCCAGCAGGCATTGATGTTACGTCATCAGGCGGCAAACTCGTCGCGAGATGTGCGTCGTCAGCTGGACAGTGGCaaa                                                                                                           >  1:1024033/1‑130 (MQ=255)
                                      gaagaaAGAGATTGGCGAGCTGGAAAACAAATTGAGCGAAACACGCGCTCGCCAGCAGGCATTGATGTTACGTCATCAGGCGGCAAACTCGTCGCGAGATGTGCGTCGTCAGCTGGACAGTGGCAAACTGGa                                                                                                      <  2:1286582/132‑1 (MQ=255)
                                      gaagaaAGAGATTGGCGAGCTGGAAAACAAATTGAGCGAAACACGCGCTCGCCAGCAGGCATTGATGTTACGTCATCAGGCGGCAAACTCGTCGCGAGATGTGCGTCGTCAGCTGGACAGTGGCAAACTGGa                                                                                                      <  2:1208632/132‑1 (MQ=255)
                                      gaagaaAGAGATTGGCGAGCTGGAAAACAAATTGAGCGAAACACGCGCTCGCCAGCAGGCATTGATGTTACGTCATCAGGCGGCAAACTCGTCGCGAGATGTGCGTCGTCAGCTGGACAGTGGCAAACTGGa                                                                                                      >  1:1208632/1‑132 (MQ=255)
                                      gaagaaAGAGATTGGCGAGCTGGAAAACAAATTGAGCGAAACACGCGCTCGCCAGCAGGCATTGATGTTACGTCATCAGGCGGCAAACTCGTCGCGAGATGTGCGTCGTCAGCTGGACAGTGGCAAACTGGa                                                                                                      >  1:1286582/1‑132 (MQ=255)
                                       aagaaagaGATTGGCGAGCTGGAAAACAAATTGAGCGAAACACGCGCTCGCCAGCAGGCATTGATGTTACGTCATCAGGCGGCAAACTCGTCGCGAGATGTGCGTCGTCAGCTGGACAGTGGCAAACTGGATGAAGCAAt                                                                                             <  1:404976/140‑1 (MQ=255)
                                        agaaagaGATTGGCGAGCTGGAAAACAAATTGAGCGAAACACGCGCTCGCCAGCAGGCATTGATGTTACGTCATCAGGCGGCAAACTCGTCGCGAGATGTGCGTCGTCAGCTGGACAGTGGCAAACTGGATGAAGCAATg                                                                                            >  2:1072547/1‑140 (MQ=255)
                                           aagaGATTGGCGAGCTGGAAAACAAATTGAGCGAAACACGCGCTCGCCAGCAGGCATTGATGTTACGTCATCAGGCGGCAAACTCGTCGCGAGATGTGCGTCGTCAGCTGGACAGTGGCAAACTGGATGAAGCAATGGCt                                                                                         <  2:296715/140‑1 (MQ=255)
                                           aagaGATTGGCGAGCTGGAAAACAAATTGAGCGAAACACGCGCTCGCCAGCAGGCATTGATGTTACGTCATCAGGCGGCAAACTCGTCGCGAGATGTGCGTCGTCAGCTGGACAGTGGCAAACTGGATGAAGCAATGGCt                                                                                         <  2:1434323/140‑1 (MQ=255)
                                           aagaGATTGGCGAGCTGGAAAACAAATTGAGCGAAACACGCGCTCGCCAGCAGGCATTGATGTTACGTCATCAGGCGGCAAACTCGTCGCGAGATGTGCGTCGTCAGCTGGACAGTGGCAAACTGGATGAAGCAATGGCt                                                                                         <  1:988648/140‑1 (MQ=255)
                                            agagATTGGCGAGCTGGAAAACAAATTGAGCGAAACACGCGCTCGCCAGCAGGCATTGATGTTACGTCATCAGGCGGCAAACTCGTCGCGAGATGTGCGTCGTCAGCTGGACAGTGGCAAACTGGATGAAGCAATGGCTc                                                                                        >  1:156578/1‑140 (MQ=255)
                                                aTTGGCGAGCTGGAAAACAAATTGAGCGAAACACGCGCTCGCCAGCAGGCATTGATGTTACGTCATCAGGCGGCAAACTCGTCGCGAGATGTGCGTCGTCAGCTGGACAGTGGCAAACTGGATGAAGCAATGGCTCGttt                                                                                    >  2:226915/1‑140 (MQ=255)
                                                    gCGAGCTGGAAAACAAATTGAGCGAAACACGCGCTCGCCAGCAGGCATTGATGTTACGTCATCAGGCGGCAAACTCGTCGCGAGATGTGCGTCGTCAGCTGGACAGTGGCAAACTGGATGAAGCAATGGCTCGTTTCGaa                                                                                >  1:75809/1‑140 (MQ=255)
                                                    gCGAGCTGGAAAACAAATTGAGCGAAACACGCGCTCGCCAGCAGGCATTGATGTTACGTCATCAGGCGGCAAACTCGTCGCGAGATGTGCGTCGTCAGCTGGACAGTGGCAAACTGGATGAAGCAATGGCTCGTTTCGaa                                                                                >  2:432891/1‑140 (MQ=255)
                                                         cTGGAAAACAAATTGAGCGAAACACGCGCTCGCCAGCAGGCATTGATGTTACGTCATCAGGCGGCAAACTCGTCGCGAGATGTGCGTCGTCAGCTGGACAGTGGCAAACTGGATGAAGCAATGGCTCGTTTCGAATCttt                                                                           >  2:1533804/1‑140 (MQ=255)
                                                          tGGAAAACAAATTGAGCGAAACACGCGCTCGCCAGCAGGCATTGATGTTACGTCATCAGGCGGCAAACTCGTCGCGAGATGTGCGTCGTCAGCTGGACAGTGGCAAACTGGATGAAGCAATGGCTCGTTTCGAATCTTTc                                                                          <  1:1265373/140‑1 (MQ=255)
                                                          tGGAAAACAAATTGAGCGAAACACGCGCTCGCCAGCAGGCATTGATGTTACGTCATCAGGCGGCAAACTCGTCGCGAGATGTGCGTCGTCAGCTGGACAGTGGCAAACTGGATGAAGCAATGGCTCGTTTCGAATCTTTc                                                                          <  2:768470/140‑1 (MQ=255)
                                                            gAAAACAAATTGAGCGAAACACGCGCTCGCCAGCAGGCATTGATGTTACGTCATCAGGCGGCAAACTCGTCGCGAGATGTGCGTCGTCAGCTGGACAGTGGCAAACTGGATGAAGCAATGGCTCGTTTCGAATCTTTCGa                                                                        <  1:131827/140‑1 (MQ=255)
                                                              aaaCAAATTGAGCGAAACACGCGCTCGCCAGCAGGCATTGATGTTACGTCATCAGGCGGCAAACTCGTCGCGAGATGTGCGTCGTCAGCTGGACAGTGGCAAACTGGATGAAGCAATGGCTCGTTTCGAATCTTTCGAAc                                                                      >  1:974967/1‑140 (MQ=255)
                                                               aaCAAATTGAGCGAAACACGCGCTCGCCAGCAGGCATTGATGTTACGTCATCAGGCGGCAAACTCGTCGCGAGATGTGCGTCGTCAGCTGGACAGTGGCAAACTGGATGAAGCAATGGCTCGTTTCGAATCTTTCGAAcg                                                                     <  2:417933/140‑1 (MQ=255)
                                                               aaCAAATTGAGCGAAACACGCGCTCGCCAGCAGGCATTGATGTTACGTCATCAGGCGGCAAACTCGTCGCGAGATGTGCGTCGTCAGCTGGACAGTGGCAAACTGGATGAAGCAATGGCTCGTTTCGAATCTTTCGAAcg                                                                     <  2:156578/140‑1 (MQ=255)
                                                               aaCAAATTGAGCGAAACACGCGCTCGCCAGCAGGCATTGATGTTACGTCATCAGGCGGCAAACTCGTCGCGAGATGTGCGTCGTCAGCTGGACAGTGGCAAACTGGATGAAGCAATGGCTCGTTTCGAATCTTTCGAAcg                                                                     <  2:1185003/140‑1 (MQ=255)
                                                               aaCAAATTGAGCGAAACACGCGCTCGCCAGCAGGCATTGATGTTACGTCATCAGGCGGCAAACTCGTCGCGAGATGTGCGTCGTCAGCTGGACAGTGGCAAACTGGATGAAGCAATGGCTCGTTTCGAATCTTTCGAAcg                                                                     <  1:1635350/140‑1 (MQ=255)
                                                               aaCAAATTGAGCGAAACACGCGCTCGCCAGCAGGCATTGATGTTACGTCATCAGGCGGCAAACTCGTCGCGAGATGTGCGTCGTCAGCTGGACAGTGGCAAACTGGATGAAGCAATGGCTCGTTTCGAATCTTTCGAAcg                                                                     <  2:1278457/140‑1 (MQ=255)
                                                               aaCAAATTGAGCGAAACACGCGCTCGCCAGCAGGCATTGATGTTACGTCATCAGGCGGCAAACTCGTCGCGAGATGTGCGTCGTCAGCTGGACAGTGGCAAACTGGATGAAGCAATGGCTCGTTTCGAATCTTTCGAAcg                                                                     <  1:1072547/140‑1 (MQ=255)
                                                               aaCAAATTGAGCGAAACACGCGCTCGCCAGCAGGCATTGATGTTACGTCATCAGGCGGCAAACTCGTCGCGAGATGTGCGTCGTCAGCTGGACAGTGGCAAACTGGATGAAGCAATGGCTCGTTTCGAATCTTTCGAAcg                                                                     <  1:1039323/140‑1 (MQ=255)
                                                               aaCAAATTGAGCGAAACACGCGCTCGCCAGCAGGCATTGATGTTACGTCATCAGGCGGCAAACTCGTCGCGAGATGTGCGTCGTCAGCTGGACAGTGGCAAACTGGATGAAGCAATGGCTCGTTTCGAATCTTTCGAAcg                                                                     <  1:1396906/140‑1 (MQ=255)
                                                                 cAAATTGAGCGAAACACGCGCTCGCCAGCAGGCATTGATGTTACGTCATCAGGCGGCAAACTCGTCGCGAGATGTGCGTCGTCAGCTGGACAGTGGCAAACTGGATGAAGCAATGGCTCGTTTCGAATCTTTCGAAcgtc                                                                   <  1:432891/140‑1 (MQ=255)
                                                                        aGCGAAACACGCGCTCGCCAGCAGGCATTGATGTTACGTCATCAGGCGGCAAACTCGTCGCGAGATGTGCGTCGTCAGCTGGACAGTGGCAAACTGGATGAAGCAATGGCTCGTTTCGAATCTTTCGAACGTCGTATTGa                                                            <  2:540267/140‑1 (MQ=255)
                                                                        aGCGAAACACGCGCTCGCCAGCAGGCATTGATGTTACGTCATCAGGCGGCAAACTCGTCGCGAGATGTGCGTCGTCAGCTGGACAGTGGCAAACTGGATGAAGCAATGGCTCGTTTCGAATCTTTCGAACGTCGTATTGa                                                            >  1:871845/1‑140 (MQ=255)
                                                                         gCGAAACACGCGCTCGCCAGCAGGCATTGATGTTACGTCATCAGGCGGCAAACTCGTCGCGAGATGTGCGTCGTCAGCTGGACAGTGGCAAACTGGATGAAGCAATGGCTCGTTTCGAATCTTTCGAACGTCGTATTGAc                                                           <  1:431530/140‑1 (MQ=255)
                                                                          cGAAACACGCGCTCGCCAGCAGGCATTGATGTTACGTCATCAGGCGGCAAACTCGTCGCGAGATGTGCGTCGTCAGCTGGACAGTGGCAAACTGGATGAAGCAATGGCTCGTTTCGAATCTTTCGAACGTCGTATTGAcc                                                          >  1:1127125/1‑140 (MQ=255)
                                                                            aaaCACGCGCTCGCCAGCAGGCATTGATGTTACGTCATCAGGCGGCAAACTCGTCGCGAGATGTGCGTCGTCAGCTGGACAGTGGCAAACTGGATGAAGCAATGGCTCGTTTCGAATCTTTCGAACGTCGTATTGACCAg                                                        <  2:1620459/140‑1 (MQ=255)
                                                                               cacGCGCTCGCCAGCAGGCATTGATGTTACGTCATCAGGCGGCAAACTCGTCGCGAGATGTGCGTCGTCAGCTGGACAGTGGCAAACTGGATGAAGCAATGGCTCGTTTCGAATCTTTCGAACGTCGTATTGACCAGATg                                                     <  1:63997/140‑1 (MQ=255)
                                                                                 cgcgCTCGCCAGCAGGCATTGATGTTACGTCATCAGGCGGCAAACTCGTCGCGAGATGTGCGTCGTCAGCTGGACAGTGGCAAACTGGATGAAGCAATGGCTCGTTTCGAATCTTTCGAACGTCGTATTGACCAGATgga                                                   <  2:958731/140‑1 (MQ=255)
                                                                                       cGCCAGCAGGCATTGATGTTACGTCATCAGGCGGCAAACTCGTCGCGAGATGTGCGTCGTCAGCTGGACAGTGGCAAACTGGATGAAGCAATGGCTCGTTTCGAATCTTTCGAACGTCGTATTGACCAGATggaagcgga                                             >  2:679958/1‑140 (MQ=255)
                                                                                       cGCCAGCAGGCATTGATGTTACGTCATCAGGCGGCAAACTCGTCGCGAGATGTGCGTCGTCAGCTGGACAGTGGCAAACTGGATGAAGCAATGGCTCGTTTCGAATCTTTCGAACGTCGTATTGACCAGATggaagcgga                                             <  1:1533804/140‑1 (MQ=255)
                                                                                             cagCCATTGATGTTACGTCATCAGGCGGCAAACTCGTCGCGAGATGTGCGTCGTCAGCTGGACAGTGGCAAACTGGATGAAGCAATGGCTCGTTTCGAATCTTTCGAACGTCGTATTGACCAGATGGAAGCGGAAGCAGa                                       <  2:500556/140‑1 (MQ=255)
                                                                                              agGCATTGATGTTACGTCATCAGGCGGCAAACTCGTCGCGAGATGTGCGTCGTCAGCTGGACAGTGGCAAACTGGATGAAGCAATGGCTCGTTTCGAATCTTTCGAACGTCGTATTGACCAGATGGAAGCGGAAGCAGaa                                      <  1:1129355/140‑1 (MQ=255)
                                                                                                gCATTGATGTTACGTCATCAGGCGGCAAACTCGTCGCGAGATGTGCGTCGTCAGCTGGACAGTGGCAAACTGGATGAAGCAATGGCTCGTTTCGAATCTTTCGAACGTCGTATTGACCAGATGGAAGCGGAAGCAGAAAg                                    >  1:1048151/1‑140 (MQ=255)
                                                                                                    tGATGTTACGTCATCAGGCGGCAAACTCGTCGCGAGATGTGCGTCGTCAGCTGGACAGTGGCAAACTGGATGAAGCAATGGCTCGTTTCGAATCTTTCGAACGTCGTATTGACCAGATGGAAGCGGAAGCAGaa                                      <  2:760961/134‑1 (MQ=255)
                                                                                                    tGATGTTACGTCATCAGGCGGCAAACTCGTCGCGAGATGTGCGTCGTCAGCTGGACAGTGGCAAACTGGATGAAGCAATGGCTCGTTTCGAATCTTTCGAACGTCGTATTGACCAGATGGAAGCGGAAGCAGaa                                      >  1:760961/1‑134 (MQ=255)
                                                                                                          tACGTCATCAGGCGGCAAACTCGTCGCGAGATGTGCGTCGTCAGCTGGACAGTGGCAAACTGGATGAAGCAATGGCTCGTTTCGAATCTTTCGAACG‑‑GTATTGACCAGATGGAAGCGGAAGCAGAAAGCCACAGCTTCgg                        <  2:974967/140‑1 (MQ=255)
                                                                                                                atcaGGCGGCAAACTCGTCGCGAGATGTGCGTCGTCAGCTGGACAGTGGCAAACTGGATGAAGCAATGGCTCGTTTCGAATCTTTCGAACGTCGTATTGACCAGATGGAAGCGGAAGCAGAAAGCCACAGCTTCGGTaaa                    >  2:1345402/1‑140 (MQ=255)
                                                                                                                atcaGGCGGCAAACTCGTCGCGAGATGTGCGTCGTCAGCTGGACAGTGGCAAACTGGATGAAGCAATGGCTCGTTTCGAATCTTTCGAACGTCGTATTGACCAGATGGAAGCGGAAGCAGAAAGCCACAGCTTCGGTaaa                    >  1:1578271/1‑140 (MQ=255)
                                                                                                                 tcaGGCGGCAAACTCGTCGCGAGATGTGCGTCGTCAGCTGGACAGTGGCAAACTGGATGAAGCAATGGCTCGTTTCGAATCTTTCGAACGTCGTATTGACCAGATGGAAGCGGAAGCAGAAAGCCACAGCTTCGGTAAAc                   >  1:754408/1‑140 (MQ=255)
                                                                                                                    ggcggcAAACTCGTCGCGAGATGTGCGTCGTCAGCTGGACAGTGGCAAACTGGATGAAGCAATGGCTCGTTTCGAATCTTTCGAACGTCGTATTGACCAGATGGAAGCGGAAGCAGAAAGCCACAGCTTCGGTAAACaaa                >  2:971804/1‑140 (MQ=255)
                                                                                                                     gcggcAAACTCGTCGCGAGATGTGCGTCGTCAGCTGGACAGTGGCAAACTGGATGAAGCAATGGCTCGTTTCGAATCTTTCGAACGTCGTATTGACCAGATGGAAGCGGAAGCAGAAAGCCACAGCTTCGGTAAACaaaa               >  1:961371/1‑140 (MQ=255)
                                                                                                                     gcggcAAACTCGTCGCGAGATGTGCGTCGCAGCTGGACAGTGGCAAACTGGATGAAGCAATGGCTCGTTTCGAATCTTTCGAACGTCGTATTGACCAGATGGAAGCGGAAGCAGAAAGCCACAGCTTCGGTAAACaaaaa              >  1:1343520/1‑140 (MQ=255)
                                                                                                                      cggcAAACTCGTCGCGAGATGTGCGTCGTCAGCTGGACAGTGGCAAACTGGATGAAGCAATGGCTCGTTTCGAATCTTTCGAACGTCGTATTGACCAGATGGAAGCGGAAGCAGAAAGCCACAGCTTCGGTAAACaaaaa              >  2:1003582/1‑140 (MQ=255)
                                                                                                                      cggcAAACTCGTCGCGAGATGTGCGTCGTCAGCTGGACAGTGGCAAACTGGATGAAGCAATGGCTCGTTTCGAATCTTTCGAACGTCGTATTGACCAGATGGAAGCGGAAGCAGAAAGCCACAGCTTCGGTAAACaaaaa              >  1:666913/1‑140 (MQ=255)
                                                                                                                         cAAACTCGTCGCGCGATGTGCGTCGTCAGCTGGACAGTGGCAAACTGGATGAAGCAATGGCTCGTTTCGAATCTTTCGAACGTCGTATTGACCAGATGGAAGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCg           >  2:270794/1‑140 (MQ=255)
                                                                                                                          aaaCTCGTCGCGAGATGTGCGTCGTCAGCTGGACAGTGGCAAACTGGATGAAGCAATGGCTCGTTTCGAATCTTTCGAACGTCGTATTGACCAGATGGAAGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGc          >  1:101431/1‑140 (MQ=255)
                                                                                                                               cgtcgCGAGATGTGCGTCGTCAGCTGGACAGTGGCAAACTGGATGAAGCAATGGCTCGTTTCGAATCTTTCGAACGTCGTATTGACCAGATGGAAGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGAc     <  1:1003582/140‑1 (MQ=255)
                                                                                                                                  cgcgagATGTGCGTCGTCAGCTGGACAGTGGCAAACTGGATGAAGCAATGGCTCGTTTCGAATCTTTCGAACGTCGTATTGACCAGATGGAAGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGAt  <  1:226915/140‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                       
CATGAAGTGACGCTGGTGGACGATACGCTGGCACGCATGAAGAAAGAGATTGGCGAGCTGGAAAACAAATTGAGCGAAACACGCGCTCGCCAGCAGGCATTGATGTTACGTCATCAGGCGGCAAACTCGTCGCGCGATGTGCGTCGTCAGCTGGACAGTGGCAAACTGGATGAAGCAATGGCTCGTTTCGAATCTTTCGAACGTCGTATTGACCAGATGGAAGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGAT  >  NC_000913/1368394‑1368663

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