Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,252,715 A→G S821P (TCC→CCC)  b1202 ← putative adhesion and penetration protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,252,7150AG83.5% 567.5 / 109.0 260S821P (TCC→CCC) b1202putative adhesion and penetration protein
Reads supporting (aligned to +/- strand):  ref base A (24/19);  new base G (126/91);  total (150/110)
Fisher's exact test for biased strand distribution p-value = 8.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.29e-01

Only 100 of 263 total aligned reads displayed.

GAGGCCATTACTCGCCTTCATAGTCATCTCATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGACAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCGTTGGCGTTAACGCCGTTGTTCTGACTGTCCAGTAC  >  NC_000913/1252623‑1252805
                                                                                            |                                                                                          
  ggCCATTACTCGCCTTCATAGTCATCTCATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGACa                                                                                          <  2:5372723/93‑1 (MQ=255)
   gCCATTACTCGCCTTCATAGTCATCTCATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAg                                                                                         <  1:4806261/93‑1 (MQ=255)
   gCCATTACTCGCCTTCATAGTCATCTCATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAg                                                                                         <  2:3099767/93‑1 (MQ=255)
   gCCATTACTCGCCTTCATAGTCATCTCATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAg                                                                                         >  2:1759594/1‑93 (MQ=255)
   gCCATTACTCGCCTTCATAGTCATCTCATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAg                                                                                         <  1:215425/93‑1 (MQ=255)
   gCCATTACTCGCCTTCATAGTCATCTCATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGACAg                                                                                         >  2:4724456/1‑93 (MQ=255)
    ccATTACTCGCCTTCATAGTCATCTCATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAgg                                                                                        >  1:4412923/1‑93 (MQ=255)
     cATTACTCGCCTTCATAGTCATCTCATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCGGGt                                                                                       >  2:4889200/1‑93 (MQ=255)
      aTTACTCGCCTTCATAGTCATCTCATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGtt                                                                                      >  1:6441935/1‑93 (MQ=255)
       ttACTCGCCTTCATAGTCATCTCATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGACAGGTTg                                                                                     >  1:3642254/1‑93 (MQ=255)
        tACTCGCCTTCATAGTCATCTCATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGa                                                                                    >  1:6909509/1‑93 (MQ=255)
        tACTCGCCTTCATAGTCATCTCATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGa                                                                                    >  1:7339803/1‑93 (MQ=255)
         aCTCGCCTTCATAGTCATCTCATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGaa                                                                                   <  2:1389113/93‑1 (MQ=255)
           tCGCCTTCATAGTCATCTCATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAAcc                                                                                 >  2:2636296/1‑93 (MQ=255)
             gCCTTCATAGTCATCTCATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACctc                                                                               <  2:7155147/93‑1 (MQ=255)
             gCCTTCATAGTCATCTCATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACctc                                                                               <  1:500054/93‑1 (MQ=255)
             gCCTTCATAGTCATCTCATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGACAGGTTGAACctc                                                                               <  1:1339971/93‑1 (MQ=255)
              ccTTCATAGTCATCTCATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACctct                                                                              >  1:659311/1‑93 (MQ=255)
                ttCATAGTCATCTCATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGc                                                                            >  2:7317448/1‑93 (MQ=255)
                ttCATAGTCATCTCATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGc                                                                            >  1:6686123/1‑93 (MQ=255)
                ttCATAGTCATCTCATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGc                                                                            >  1:3473593/1‑93 (MQ=255)
                ttCATAGTCATCTCATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGc                                                                            >  2:1676202/1‑93 (MQ=255)
                 tCATAGTCATCTCATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGcc                                                                           <  1:6122145/93‑1 (MQ=255)
                 tCATAGTCATCTCATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGcc                                                                           >  1:2590230/1‑93 (MQ=255)
                 tCATAGTCATCTCATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGACAGGTTGAACCTCTGcc                                                                           >  2:1105654/1‑93 (MQ=255)
                     aGTCATCTCATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGc                                                                       <  2:4295458/93‑1 (MQ=255)
                      gTCATCTCTTTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAAAACCCATAACCAGAAGGGGACAGGTTGAACCTCTCCCCGCCt                                                                      >  1:2376183/1‑93 (MQ=255)
                      gTCATCTCATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCt                                                                      >  1:6182553/1‑93 (MQ=255)
                      gTCATCTCATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCt                                                                      >  1:4568436/1‑93 (MQ=255)
                      gTCATCTCATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGACAGGTTGAACCTCTGCCCGGCt                                                                      >  1:2658771/1‑93 (MQ=255)
                       tCATCTCATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCtt                                                                     >  1:7446313/1‑93 (MQ=255)
                        cATCTCATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTc                                                                    >  2:618038/1‑93 (MQ=255)
                        cATCTCATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTc                                                                    >  1:6046275/1‑93 (MQ=255)
                         aTCTCATTCTGGTGGCTGTATGTAAGCTGGGTTTGTGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTcc                                                                   >  1:6994236/1‑93 (MQ=255)
                         aTCTCATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTcc                                                                   <  2:5136717/93‑1 (MQ=255)
                          tctcATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCa                                                                  >  1:5610484/1‑93 (MQ=255)
                          tctcATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCa                                                                  <  2:384766/93‑1 (MQ=255)
                            tcATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAgg                                                                <  1:4931456/93‑1 (MQ=255)
                            tcATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAgg                                                                >  2:4018299/1‑93 (MQ=255)
                            tcATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAgg                                                                >  2:6596166/1‑93 (MQ=255)
                            tcATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAgg                                                                <  2:2172324/93‑1 (MQ=255)
                            tcATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGACAGGTTGAACCTCTGCCCGGCTTCCAgg                                                                <  1:4256875/93‑1 (MQ=255)
                              aTTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAGCCTCTGCCCGGCTTCCAGGga                                                              <  2:7249733/93‑1 (MQ=255)
                              aTTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGga                                                              <  1:5243719/93‑1 (MQ=255)
                               ttCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGgag                                                             >  1:5471648/1‑93 (MQ=255)
                               ttCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGACAGGTTGAACCTCTGCCCGGCTTCCAGGgag                                                             >  1:1473561/1‑93 (MQ=255)
                                tCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGgaga                                                            <  1:1103393/93‑1 (MQ=255)
                                tCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGgaga                                                            <  2:3418559/93‑1 (MQ=255)
                                tCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGgaga                                                            <  2:2895585/93‑1 (MQ=255)
                                tCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGACAGGTTGAACCTCTGCCCGGCTTCCAGGgaga                                                            <  1:3375438/93‑1 (MQ=255)
                                 cTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGAt                                                           <  1:885934/93‑1 (MQ=255)
                                  tggtggCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATg                                                          >  1:2683803/1‑93 (MQ=255)
                                  tggtggCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGACAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATg                                                          >  1:5647702/1‑93 (MQ=255)
                                     tggCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTc                                                       <  2:6305706/93‑1 (MQ=255)
                                       gCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTcat                                                     >  2:5708924/1‑93 (MQ=255)
                                         tgtatgtaAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTcattc                                                   >  1:71709/1‑93 (MQ=255)
                                          gtatgtaAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTcattcc                                                  >  2:3558226/1‑93 (MQ=255)
                                             tgtaAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTcattccatt                                               >  1:1517747/1‑93 (MQ=255)
                                             tgtaAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGACAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTcattccatt                                               <  2:5909790/93‑1 (MQ=255)
                                             tgtaAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGACAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTcattccatt                                               <  1:7501201/93‑1 (MQ=255)
                                              gtaAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTcattccattc                                              >  2:6655402/1‑93 (MQ=255)
                                              gtaAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTcattccattc                                              <  2:5237659/93‑1 (MQ=255)
                                                aaGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGc                                            >  2:199939/1‑93 (MQ=255)
                                                 aGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCa                                           <  2:3212551/93‑1 (MQ=255)
                                                 aGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCa                                           >  2:5942253/1‑93 (MQ=255)
                                                   cTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGt                                         <  1:3278539/93‑1 (MQ=255)
                                                    tGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGACAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTg                                        >  1:1605640/1‑93 (MQ=255)
                                                     gggTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGACAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGc                                       >  1:4498006/1‑93 (MQ=255)
                                                      ggTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGcc                                      >  1:3854336/1‑93 (MQ=255)
                                                      ggTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGACAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGcc                                      <  1:6530272/93‑1 (MQ=255)
                                                       gTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCg                                     <  1:1833816/93‑1 (MQ=255)
                                                       gTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCg                                     <  2:1170444/93‑1 (MQ=255)
                                                       gTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCg                                     <  2:7381685/93‑1 (MQ=255)
                                                        tttGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCGt                                    >  1:193023/1‑93 (MQ=255)
                                                          tGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCGTTg                                  >  2:5148142/1‑93 (MQ=255)
                                                          tGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCGTTg                                  <  1:3095177/93‑1 (MQ=255)
                                                           ggggCTCTATATAGACCCCATACCAGGTAGGGGGCAGGTTGACCCTCTGCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCGTTGGc                                <  2:7109543/91‑1 (MQ=255)
                                                              gCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCGTTGGCGt                              >  2:814569/1‑93 (MQ=255)
                                                              gCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCGTTGGCGt                              >  2:1886649/1‑93 (MQ=255)
                                                               ctctATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCGTTGGCGtt                             >  1:346451/1‑93 (MQ=255)
                                                                tctATATAGAACCCATAACCAGTAGGGGACAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCGTTGGCGTTa                            >  1:293326/1‑93 (MQ=255)
                                                                  tatataGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCGTTGGCGTTAAc                          <  2:224081/93‑1 (MQ=255)
                                                                  tatataGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCGTTGGCGTTAAc                          >  1:6314103/1‑93 (MQ=255)
                                                                   atataGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCGTTGGCGTTAACg                         <  1:1893451/93‑1 (MQ=255)
                                                                         aaCCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCGTTGGCGTTAACGCCgttg                   >  2:225096/1‑93 (MQ=255)
                                                                          aCCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCGTTGGCGTTAACGCCgttgt                  >  1:1394964/1‑93 (MQ=255)
                                                                            ccATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCGTTGGCGTTAACGCCGTTGTTc                <  1:4048506/93‑1 (MQ=255)
                                                                              aTAACCAGTAGGGGACAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCGTTGGCGTTAACGCCGTTGTTCTg              <  1:2294148/93‑1 (MQ=255)
                                                                               tAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCGTTGGCGTTAACGCCGTTGTTCTGa             >  2:266348/1‑93 (MQ=255)
                                                                               tAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCGTTGGCGTTAACGCCGTTGTTCTGa             >  1:3975629/1‑93 (MQ=255)
                                                                                 aCCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCGTTGGCGTTAACGCCGTTGTTCTGACt           <  1:2829533/93‑1 (MQ=255)
                                                                                    aGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCGTTGGCGTTAACGCCGTTGTTCTGACTgnn        >  2:2270966/1‑91 (MQ=255)
                                                                                    aGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCGTTGGCGTTAACGCCGTTGTTCTGACTGTc        <  1:5768988/93‑1 (MQ=255)
                                                                                    aGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCGTTGGCGTTAACGCCGTTGTTCTGACTGTc        <  2:4553859/93‑1 (MQ=255)
                                                                                     gTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCGTTGGCGTTAACGCCGTTGTTCTGACTGTcc       <  2:6413212/93‑1 (MQ=255)
                                                                                     gTAGGGGACAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCGTTGGCGTTAACGCCGTTGTTCTGACTGTcc       <  2:4201337/93‑1 (MQ=255)
                                                                                      tAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCGTTGGCGTTAACGCCGTTGTTCTGACTGTCCa      <  2:6820543/93‑1 (MQ=255)
                                                                                      tAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCGTTGGCGTTAACGCCGTTGTTCTGACTGTCCa      >  2:2558642/1‑93 (MQ=255)
                                                                                        ggggACAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCGTTGGCGTTAACGCCGTTGTTCTGACTGTCCAGt    >  2:1475692/1‑93 (MQ=255)
                                                                                         ggggCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCGTTGGCGTTAACGCCGTTGTTCTGACTGTCCAGTa   <  1:1245580/93‑1 (MQ=255)
                                                                                            |                                                                                          
GAGGCCATTACTCGCCTTCATAGTCATCTCATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGACAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCGTTGGCGTTAACGCCGTTGTTCTGACTGTCCAGTAC  >  NC_000913/1252623‑1252805

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